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ABSTRACT 

We present maps of the CO-to-H2 conversion factor (aco) and dust-to-gas ratio (DGR) in 26 nearby, 
star-forming galaxies with ^kiloparsec spatial resolution. We have simultaneously solved for aco and 
DGR by assuming that the DGR is approximately constant on kpc scales. With this assumption, we 
can combine maps of dust mass surface density, CO integrated intensity and HI column density to 
solve for both aco and DGR with no assumptions about their value or dependence on metallicity 
or other parameters. Such a study has just become possible with the availability of high resolution 
far-IR maps from the Herschel key program KINGFISH, 12 CO J=(2— 1) maps from the IRAM 30m 
large program HERACLES and HI 21-cm line maps from THINGS. We use a fixed ratio between the 
(2—1) and (1—0) lines to present our aco results on the more typically used 12 CO J=(l— 0) scale 
and show using literature measurements that variations in the line ratio do not effect our results. In 
total, we derive 782 individual solutions for aco and DGR. On average, aco =3.1 M pc -2 (K km 
s -ij-i £ Qr Qur sam pi e w ith a standard deviation of 0.3 dex. Within galaxies we observe a generally 
flat profile of aco as a function of galactocentric radius. However, most galaxies exhibit a lower aco 
in the central kiloparsec — a factor of ~2 below the galaxy mean, on average. In some cases, the central 
aco value can be factors of 5 to 10 below the standard MW value of aco,MW = 4.4 M pc -2 (K 
km s -1 ) -1 . While for aco we find only weak correlations with metallicity, DGR is well-correlated 
with metallicity, with an approximately linear slope. Finally, we present several recommendations for 
choosing an appropriate aco fc> r studies of nearby galaxies. 

Keywords: dust, extinction — infrared: ISM — ISM: ISM — molecules: Galaxies — ISM 
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1. INTRODUCTION 

The H2 molecule is difficult to observe in the preva- 
lent interstellar medium (ISM) conditions of a normal 
star-forming galaxy. Since it is the primary constituent 
of molecular gas, inferring the mass of H2 is crucial for 
studying this phase of the ISM and the star formation 
that occurs within it. The widespread practice is to use 
the second most abundant molecule, carbon monoxide 
( 12 CO), as a tracer and convert measured CO integrated 
intensities into H2 column densities using a "CO-to-Hy 
conversion factor XcoQ 



face density with information on the DGR, i.e. using the 
following equation: 



Nh 2 = Xcolco- 



(i) 



In mass surface density units this equation can be rewrit- 
ten as 

Sh 2 = acolco (2) 

where Eh 2 is the total mass surface density of molecular 
gas (including a correction for the second most abundant 
element, He). A variety of observations have shown that 
ftco~ 4.4 M pc -2 (K km s " 1 ) -1 is characteristic of the 
local area of the Milky Way (Solomon et al.|1987 Strong 
fc Mattox|1996| |Abdo et al.|2UlUp . Despite uncertainties 
in the physics behind the conversion factor, the observ- 
ability of CO ensures that it will remain a widely- used 
molecular gas tracer, particularly at high redshift. For 
this reason, measurements of aco in nearby galaxies, 
where we can relate its variations to local environmental 
conditions, are critical. 

In order to measure the conversion factor, one must 
measure Eh 2 (or equivalently, E gas — Ehi) independently 
of CO and then compare it to observed CO integrated 
intensities. This has been done using a variety of tech- 
niques, including: measuring total gas masses from 7- 
ray emission plus a mode l for the cosmic ray distribution 



(Strong fc Mattox 1996 Abdo et al. 2010), using the 



observed velocity dispersion and size of the molecular 
cloud to obtain virial masses (ISolomon et al.||1987l IWiT 
son|1995||Bolatto et al.|2008HDonoyan Meyer eTaT"| 20 12 ; 
Wei et al.| |2012| |Gratier et al.||2012p , and modeling mul- 



tiple molecular gas lin es with varying op tical d epths and 
critical densities (e.g. |Weifi et al.||200l| |Israel||2C)09b|a| >. 
In general, few of these techniques are effective tor con- 
straining aco m galaxies outside the Local Group due 
to the difficulty of obtaining the necessary observations 
(e.g. 7-ray maps or CO mapping at < lOOpc resolution) 
or doubts about fundamental assumptions (e.g. that CO 
traces the full extent of molecular gas in the clouds; that 
the clouds lack contributions to virial balance from mag- 
netic or pressure forces; that simple radiative transfer 
models can reproduce molecular gas excitation on kpc- 
scales). 

Another potential technique to measure Eh 2 is to use 
dust as a tracer of total gas column. Assuming that dust 
and gas are well mixed, the dust-to-gas ratio (DGR) is 
not a function of atomic/molecular phase and the frac- 
tion of mass in ionized gas is negligible, the observed dust 
mass surface density can be converted to a gas mass sur- 



1 We refer to the CO-to-H2 conversion factor in mass units 
throughout this paper. Including a factor of 1.36 for helium, 
X C o= 2 x 10 20 cm -2 (K km s -1 ) -1 corresponds to aco= 4.35 
Mq pc -2 (K km s -1 ) -1 . aco can t> e converted to Xco units by 
multiplying by a factor of 4.6 x 10 19 . 



DGR 



-•HI 



Eh 2 = Shi + acolco- 



(3) 



Here Ed is the mass surface density of dust and Ehi 
and Eh 2 are the mass surface densities of atomic and 
molecular ga^J respectively, and DGR is the dust-to-gas 
mass ratio. By measuring Ed and Ehi and assuming 
the DGR (or simultaneously measuring it, as we will dis- 
cuss shortly), the molecular gas mass can be determined. 
While still subject to its own systematic uncertainties, 
this technique relies on a different set of assumptions 
than those mentioned previously. 

Studies of DGR or aco have typically fixed one of the 
parameters in order to determine the other, so it is diffi- 
cult to avoid circularity when using a fixed DGR to solve 
for aco m Equation [3] Alternatively, with sufficiently 
high spatial resolution, DGR can be determined on sight- 
lines free of molecular gas and extrapolated to regions 
where CO is detected. For very nearby galaxies like the 
Magel l anic C loud s, such a technique has been used by 
Israeli fll997| ) and |Leroy et al] fl2009b| ), who found that 
aco determinations from virial masses can be strongly 
biased by envelopes of "CO-dark" molecular gas in low- 
metallicity systems. In more distant galaxies, we gen- 
erally cannot isolate purely atomic lines of sight at the 
resolution of typical HI and CO observation s. Instead, it 
is po ssible to use a technique developed by |Leroy et al. 
(2011) to simultaneously measure aco and DGR using 
the assumption that DGR should be constant over a re- 
gion of a galaxy. Since we can now achieve ~kpc spatial 
resolution in the far-IR with Herschel and have sensitive, 
high-resolution CO and HI maps, it is possible to extend 
this technique, which has thus far only been applied to 
the Local Group, to more distant galaxies. 

The idea behind this technique (explained in detail in 
Section |3| is that spatially resolved measurements of Ed, 
Ehi ana Ico allow one to solve for aco ana " DGR in 
a region when it is smaller than the typical scale over 
which DGR varies (i.e. it is well represented by one DGR 
value) and it covers a range of CO/HI ratios. Having 
chosen an appropriately sized region, we can determine 
these two constants, defining the best fit as that which 
produces the most uniform DGR for the multiple lines- 
of-sight included in the region. This technique assumes 
no prior value of DGR or aco — it makes use of the (by 
definition) linear dependence of Eh 2 on aco to identify 
the solution. The technique is therefore only applica- 
ble in regions where CO is detected. Constraining aco 
in more ext reme conditions woul d require different tech- 



2012 for more discussion). 



niques (see Schruba et al. 

One useful aspect of solving simultaneously for aco 
and DGR is that the absolute normalization of the dust 
tracer is irrelevant for the determination of aco as l° n g 
as it is linear with the "true" Ed- Since the DGR is cal- 
ibrated from the map itself, any constant term will show 
up in both DGR and Ed and makes no impact on the as- 
sessment of Eh 2 



For example, Leroy et al 



(2011) showed 



2 In converting from column densities to mass surface densities, 
we account for helium with a factor of 1.36. For Eh 2 , this factor is 
included in the aco term. We apply the helium correction to Ehi 
as well. 



aco & DGR in Nearby Galaxies 



3 



that the dust optical depth at 160 /im, determined using 
an assumed power-law dependence on frequency and an 
estimate of the dust temperature from the 70/160 ratio, 
works comparably we ll as the dust mass surface density. 
|Dobashi et al.| ([2008 ) performed a similar study in the 
Large Magellanic Cloud using Ay mapping to trace dust 
mass surface density. Since we are also interested in the 
value of DGR itself, we use Ed as our tracer throughout 
this paper. We note that uncertainties in the absolute 
value of Ed, as long as they do not introduce a non- 
linearity within the region in question, will not affect the 
determination of aco- 

For this study, we make use of CO J=(2 — 1) observa- 
tions from the large pro gram HERACLES on the IRAM 
30m ( Leroy et al.|2009a| ). It is important to note that we 
are therefore determining aco appropriate for that CO 
line. Since most studies quote ftco f° r the CO J=(l — 0) 
line, throughout the paper we use a line ratio (R21) to 
convert our measurements to the (1—0) scale. Systematic 
variations in R21 will result in errors in the (1 — 0) con- 
version factor, but the (2 — 1) conversion factor will not 
be affected since it is what we are directly deriving. Al- 
though it is convenient to discuss ftco m its typical (1—0) 
incarnation, we note that the (2 — 1) conversion factor it- 
self will be important for future studies with ALMA. For 
nearby galaxies, at a given angular resolution, mapping 
in the (2 — 1) line is more efficient than in the (1 — 0) line. 
For high redshift galaxies, the (1 — 0) line may be shifted 
out of ALMA's frequency coverage. Thus, aco for the 
(2 — 1) line will be useful regardless of its relationship to 

(i-o). 

This paper is laid out as follows: Section [2] presents the 
details of the resolved dust and gas observations we use. 
In Section [3] we describe the technique to simultaneously 
measure DGR and aco and discuss how the procedure is 
optimized to deal with more distant galaxies than those 
in the Local Group (more details on the technique can be 
found in the Appendix). We present the results of per- 
forming the solution on 26 nearby galaxies in Section [4] 
and discuss their implications for our understanding of 
how DGR and aco vary with metallicity and other ISM 
properties in Section [6] 

2. OBSERVATIONS 

We make use of observations of dust and gas from a 
series of surveys of nearby galaxies built upon SINGS 



(the "Spitzer In frared Nearby Galaxies Survey"; Kenni- 
|cutt et al.|2003| >. This includ es HI from "The HI Nearby 



Galaxies Survey" (THINGS; | Walter et al.||2008[) , 12 CO 
J=(2-l) from the "HERA CU L ine Emission Survey" 



(HERACLES; |Leroy et al.1|2009aj Leroy et. al, in prep) 
and far-infrared dust emission from "Key Insights into 
Nearby Gala xies: A Far-Infrared Su rvey with Herschel" 
(KINGFISH; |Kennicutt et al ||2011| . In addition, sev- 
eral galaxies that were not included in THINGS have HI 
observations from either archival or new observations. 
The sample of galaxies in common between these sur- 
veys and for which we h ave detections of CO emission 
(see |Schruba et aL|[2012j for details on the HERACLES 
non-detections) consists of 26 targets listed in Table HI 

The SINGS and KINGFISH surveys targeted galaxies 
with a variety of morphologies, located within 30 Mpc. 
Due to the requirement of a CO detection for our work, 
all but one of the viable targets are spiral galaxies, the 



exception being NGC 3077 which is a starbursting dwarf. 
In Table [I] we list the positions, distances, orientation 
parameters and the B-band isophotal radii at 25 mag 
arcsec -2 (R25) for our targets. Throughout the text we 
define T25 = r/ii^s, where r is the galactocentric radius 
in units of arcminutes. 



2.1. Dust Mass Surface Density 

We use dust mass surface density maps derived from 
pixel-by-pixel modeling of the infrared (IR) spectral en- 
ergy distribution (SED) o bserved by Spitzer a nd Herschel 
with models developed by |Draine fc Li| ( |2007| ). A detailed 
description of t he modeling lor NGC^ 6251md NGC 6946 
is presented in Aniano et al. (2012) and the full sam- 
ple results are presented in Aniano et al (in prep). The 
dust models include a description of the dust properties 
(size distribution, composition and optical properties of 
the grains) with a variable fraction of dust in the form 
of polycyclic aromatic hydrocarbons (PAHs). The dust 
is illuminated by a radiation field distribution wherein 
a fraction of the dust is heated by a minimum radia- 
tion field U m in while the rest is heated by a power-law 
distribution of radiation fields extending up to U =10 7 
Ummp, w here Ummp is the sol ar neighborhood radiation 
field from |Mathis et aL (1983). The fraction of the dust 
mass heated by radiation fields where U>10 2 Ummp is 
typically very small, so the dust mass surface density is 
not very sensitive to the exact value of the upper U limit. 

The resolution of the dust mass surface density map 
is equivalent to that of the lowest resolution IR map 
included in the modeling. In order to preserve spatial 
resolution while still covering the peak of the dust SED, 
we use the dust modeling at a resolution matched to the 
SPIRE 350 /im map (FWHM - 2b"). This limiting res- 
olution allows us to include the following maps in the 
dust modeling: IRAC 3.6, 4.5, 5.8 and 8.0 /im; MIPS 24 
and 70 /im; PACS 70, 100 and 160 /im; and SPIRE 250 
and 350 /im. In theory, we could perform this analysis 
at even higher resolution using the SPIRE 250 resolu- 
tion maps, which include all IRAC bands; MIPS 24 /im; 
all PACS ba nds; and SPIRE 250 /im. However, |Aniano| 
et al. (2012) found that maps where the limiting reso- 
lution exceeds MIPS 70 /im are less reliable due to the 
comparatively low surface brightness sensitivity of the 
PACS observations. At SPIRE 250 /im resolution, they 
find systematic errors of up to ~ 30 — 40% in the dust 
mass (compared to their best estimate, which includes all 
IRAC, MIPS, PACS and SPIRE bands). However, when 
both SPIRE 250 and 350 observations are included, the 
systematic errors in the dust mass are ~ 10% or less. We 
proceed by using the SPIRE 350 resolution dust model- 
ing results. 

Aside from the dust mass surface density, the |Aniano| 
et al. (2012) modeling also constrains a number of other 



quantities that we make use of later in interpreting the 
results. These include U, the average radiation field heat- 
ing the dust; U m i n the minimum radiation field; /pdr 
the fraction of the dust luminosity that comes from dust 
heated by U > 100 Ummp; and <?pah the fraction of the 
dust mass accounted for by PAHs with fewer than 10 3 
carbon atoms. 
Statistical uncertainties on Hp and the other derived 



parameters were measured by Aniano et al. (2012) with 
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Table 1 

Galaxy Sample 



Galaxy 


R.A. 


Dec. 


Distance 


i 


P.A. 


Morphology 


#25 
(') 


Solution Pixel a 




(J2000) 


(J2000) 


(Mpc) 


(°) 


(°) 




Radius (kpc) 


NGC 0337 


00 h 59 m 50!l 


-07°34 / 41 / /0 


19.3 


51 


90 


SBd 


1.5 


3.5 


NGC 0628 


01 h 36 m 41!8 


+15°47 / 00 / /0 


7.2 


7 


20 


SAc 


4.9 


1.3 


NGC 0925 


02 h 27 m 16!5 


+33°34 / 43 / /5 


9.1 


66 


287 


SABd 


5.4 


1.7 


NGC 2841 


09 h 22 m 02!6 


+50°58 / 35 / /2 


14.1 


74 


153 


SAb 


3.5 


2.6 


NGC 2976 


09 h 47 m 15!3 


+67°55 / 00 / /0 


3.6 


65 


335 


SAc 


3.6 


0.6 


NGC 3077 


10 h 03 m 19!l 


+68°44 / 02 / /0 


3.8 


46 


45 


IOpec 


2.7 


0.7 


NGC 3184* 


10 h 18 m 17!0 


+41°25'28 / /0 


11.7 


16 


179 


SABcd 


3.7 


2.1 


NGC 3198 


10 h 19 m 55!0 


+45°32 / 58 / /9 


14.1 


72 


215 


SBc 


3.2 


2.6 


NGC 3351* 


10 h 43 m 57!7 


+ll°42 / 14 / /0 


9.3 


41 


192 


SBb 


3.6 


1.7 


NGC 3521* 


ll h 05 m 48!6 


-00°02 / 09 / /2 


11.2 


73 


340 


SABbc 


4.2 


2.0 


NGC 3627* 


ll h 20 m 15!0 


+12°59 / 29 / /6 


9.4 


62 


173 


SABb 


5.1 


1.7 


NGC 3938 


ll h 52 m 49!4 


+44°07 / 14 / /9 


17.9 


14 


195 


SAc 


1.8 


3.3 


NGC 4236 


12 h 16 m 42!l 


+69°27 / 45 / /0 


4.5 


75 


162 


SBdm 


12.0 


0.8 


NGC 4254* 


12 h 18 m 49!6 


+14°24 / 59 / /0 


14.4 


32 


55 


SAc 


2.5 


2.6 


NGC 4321* 


12 h 22 m 54!9 


+15°49 / 21 / /0 


14.3 


30 


153 


SABbc 


3.0 


2.6 


NGC 4536* 


12 h 34 m 27!l 


+02°ll / 16 / /0 


14.5 


59 


299 


SABbc 


3.5 


2.6 


NGC 4569 


12 h 36 m 49!8 


+13°09 / 46 / /0 


9.9 


66 


23 


SABab 


4.6 


1.8 


NGC 4625 


12 h 41 m 52!7 


+41°16 / 26 / /0 


9.3 


47 


330 


SABmp 


0.7 


1.7 


NGC 4631 


12 h 42 m 08?0 


+32°32 / 29 / /0 


7.6 


85 


86 


SBd 


7.3 


1.4 


NGC 4725 


12 h 50 m 26!6 


+25°30 / 03 / /0 


11.9 


54 


36 


SABab 


4.9 


2.2 


NGC 4736* 


12 h 50 m 53!0 


+41°07 / 13 / /2 


4.7 


41 


296 


SAab 


3.9 


0.8 


NGC 5055* 


13 h 15 m 49!2 


+42°01 / 45 / /3 


7.9 


59 


102 


SAbc 


5.9 


1.4 


NGC 5457 b * 


14 h 03 m 12!6 


+54°20 / 57 / /0 


6.7 


18 


39 


SABcd 


12.0 


1.2 


NGC 5713 


14 h 40 m ll!5 


-00°17 / 21 / /0 


21.4 


48 


11 


SABbcp 


1.2 


3.9 


NGC 6946* 


20 h 34 m 52!2 


+60°09 / 14 / /4 


6.8 


33 


243 


SABcd 


5.7 


1.2 


NGC 7331 


22 h 37 m 04!0 


+34°24 / 56 / /5 


14.5 


76 


168 


SAb 


4.6 


2.6 



Note. 



from [Walter et al.| ( |2008| ) where possible, LEDA and NED databases otherwise. 

a Solution pixels are defined in Section [3] They are the regions in which we solve for cnco and DGR. 

£ M101 _ . 

CO J= (1 — 0) maps available from Nobeyama survey of nearby galaxies (Kuno et al. 2007). 



a Monte Carlo approach. In most of the regions we con- 
sider the statistical uncertainties on the dust mass sur- 
face densities are small, but the S/N of the dust mass 
maps is generally a function of radius and the outskirts 
can have S/N ~ 5 in some cases. Our error estimates 
take these uncertainties into account and are described 
in Section 



Table 2 

HI Observation Summary 



Galaxy 



Sou 



FWHM Beam Properties 



in Section 



3.3 



3.5 



Potential systematic effects are discussed 



2.2. HI Surface Density 

To trace the atomic gas surface density in our targets, 
we use a combi nation of NRA OE] VLA HI observations 
from THINGS ( [Walter et aL]|2008| ) and supplementary 
HI, both new and archival, described in Leroy et al. (in 
prep.). The source and angular resolution of the HI maps 
for our targets are listed in Table [2] 

The HI maps are converted from integr ated intensi- 
ties to column densities as in Equation 5 of W alter et al.| 
(2008). We convolve each map with an elliptical Gaus- 
sian kernel determined by its individual beam proper- 
ties to produce a circular Gaussian point spread func- 
tion (PSF). We then use kernels created following the 
procedures in |Aniano et al.|(|2011[ ) to convolve the circu- 
lar Gaussian to match the PSF at SPIRE 350 /im. For 
the HI we assume the uncertainties to be the larger of 

3 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 







Major (") 


Minor (") 


P.A. ( u ) 


NGC 0337 


Archival 


20.2 


13.0 


160.3 


NGC 0628 


THINGS 


11.9 


9.3 


-70.3 


NGC 0925 


THINGS 


5.9 


5.7 


30.6 


NGC 2841 


THINGS 


11.1 


9.4 


-12.3 


NGC 2976 


THINGS 


7.4 


6.4 


71.8 


NGC 3077 


THINGS 


14.3 


13.2 


60.5 


NGC 3184 


THINGS 


7.5 


6.9 


85.4 


NGC 3198 


THINGS 


11.4 


9.4 


-80.4 


NGC 3351 


THINGS 


9.9 


7.1 


24.1 


NGC 3521 


THINGS 


14.1 


11.2 


-61.7 


NGC 3938 


Archival 


18.5 


18.2 


48.6 


NGC 3627 


THINGS 


10.0 


8.9 


-60.9 


NGC 4236 


New 


16.7 


13.9 


69.6 


NGC 4254 


Archival 


16.9 


16.2 


54.4 


NGC 4321 


Archival 


14.7 


14.1 


163.4 


NGC 4536 


Archival, New 


14.7 


13.8 


-11.4 


NGC 4569 


Archival 


14.2 


13.9 


32.9 


NGC 4625 


Archival 


13.0 


12.5 


-29.2 


NGC 4631 


Archival 


14.9 


13.3 


178.1 


NGC 4725 


Archival, New 


18.6 


17.0 


-20.9 


NGC 4736 


THINGS 


10.2 


9.1 


-23.0 


NGC 5055 


THINGS 


10.1 


8.7 


-40.0 


NGC 5457 


THINGS 


10.8 


10.2 


-67.1 


NGC 5713 


Archival 


15.5 


14.9 


121.9 


NGC 6946 


THINGS 


6.0 


5.6 


6.6 


NGC 7331 


THINGS 


6.1 


5.6 


34.3 



Note. — Galaxies with new HI data were observed in VLA project 
AL735 (P.I. A. Leroy). 
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either 0.5 M pc -2 or 10% of the measured column den- 
sity. Potential systematic unce rtain ties from HI opacity 
effects are discussed in Section 13^41 



2.3. CO Integrated Intensity 

To trace the molecular gas distribution in our tar- 
gets we use CO J=(2-l) ma pping from the HERACLES 



survey QLeroy et al.||2009a| ). Integrated intensity maps 
were generated from the CO spectral cubes by integrat- 
ing the spectra over a range in velocity around either 
1) the detected CO line in that spectrum or 2) the ex- 
pected CO velocity predicted from the HI velocity (since 
HI is detected at high S/N in almost all relevant pix- 
els). We propagate uncertainties through these masking 
steps, creating in the end an integrated CO line map and 
an uncertainty map. The HERACLES maps have PSFs 
well approximated by a circular Gaussian with FWHM 
of 13 A" . We convolve the maps with kernels construc ted 
using the techniques described by |Aniano et al. (2011 ) to 
match the resolution of the SPIRE 350 maps. We have 
tested the effect of error beams (i.e. extended wings of 
the PSF or stray light pick-up of the IRAM 30-m) on 
the HERACLES maps and find that the effect is less 
than ~ 5% for all galaxies and closer to 1% for most. 

Our measurements of aco directly determine the con- 
version factor appropriate for CO J=(2 — 1). However, to 
compare with the standard CO J=(l — 0) factor, which 
is more frequently used, we convert to (1 — 0) using a 
fixed value of the line ratio i?2i = (2 — 1)/(1 — 0) = 0.7. 
Due to revised telescope efficiencie s, the R<2\ we use dif - 
fers slightly from what was used in |Leroy et al.j ( |2009a[ ) . 
The R21 we assume was found to be an appropriate aver- 
age for the HERACLES sample (E. Rosolowsky, private 
communication). We discuss the effects of assuming a 
fixed R21 on o ur results for the (1 — 0) conversion factor 
in Section |3.4| To apply the aco we report in the fol- 
lowing sections to CO J=(2 — 1) observations, its value 
should be divided by 0.7. 

2.4. Ancillary Datasets 
2.4.1. Metallicity 

In the analysis presented in Section [4] we study the 
variations of aco and DGR as a function of metallicity. 
Wherever possibl e, we make use of the metallicity mea- 
surements from |Moustakas et al. (2010 hereafter M10) 
who derived characteristic metallicities as well as radial 
gradients in oxygen abundance for the SINGS sample. 
M10 present results using two dif ferent calibrations for 
the str ong line ab undances — from [ Kobulnic ky fc Ke wley 
( |2QQ4| KK04) and |Pilyugin Thuan| ( |2005| PT05). Both 
calibrations are considered in the following analysis. 

Our preferred metallicity measurement is a radial gra- 
dient from HII region metallicities (M10, Table 8). For 
several galaxies no radial gradient measurement is avail- 
able, so we use a fixed metallicity for the entire galaxy 
equal to the "characteristic metallicity" (M10, Table 9). 
In the case of NGC 4236 and 4569, the only metallicity 
measurements available are from the B-band luminosity- 
metallicity (L — Z) relationship and we use those values 
from M10 with no gradient. Finally, two of our galax- 
ies are not in the M10 sample. For NGC 3077 we use a 
metallicit y of 12+log(Q/H)= 8 .9 in KK04 with no gradi- 
ent, from Calzetti et al.l (119941). To obtain a PT05 mea- 



surement for NGC 3077, we subtract 0.6 dex, the aver- 
age offset between the two calibrations found by M10. 
For NGC 5457 (a.k.a. M101), a galaxy with a well- 
known radia l metall i city g radient, we use the measure- 
ments from Bresolin (2007). These metallicities are from 
direct methods, so are not on the same scale as either 
PT05 or KK04. The metallicities and gradients we use 
are listed in Table 03 These are converted to match the 
r 2 5 we adopt in this work, which can be slightly different 
from that adopted by M10. In order to compare with 
the local MW, we use the metallicity of the Orion Neb- 
ula HII region in the PT05 and KK04 calibrations from 
work by Croxall et al. (in prep), i.e. 12+log(0/H)=8.5 
for PT05 and 8.8 for KK04. 

2.4.2. Star Formation Rate and Stellar Mass Surface 
Density Maps 

The star formation rate surface densities (Ssfr) are 
calculated from Ha and 24 /im maps using the Ha maps 
and procedure described in |Leroy et al. ( |2012[ ). The Ha 
maps have been convolved to match the SPIRE 350 jam 
PSF assuming an initial - 1-2" FWHM Gaussian PSFs, 
although the large difference between the initial and fi- 
nal PSF makes this choice mostly irrelevant. We use 
24 /am maps from SINGS, processed (backgrou nd sub- 
tracted, con volved and aligned) as des cribed in Aniano 
|et al.[ fl2012| ) . The |Aniano et at] ( |2012D modeling results 
described in Section \l. II are also used to remove a non- 
star-formation relate d cirrus componen t from the 24 /xm 
map as described in Leroy et al. ( |2012[ ). 

We calculate the stellar mass surface density (£*) from 
th e IRAC 3.6 / i m obs ervations from SINGS, as described 



Leroy et al. ([2008). This provides only a rough tracer 



of stellar mass surface density, since we do not take into 
account c orrections for various contaminants in the 3.6 
/xm band flZibetti fc Groves||2011| |Meidt et al.]|2012 ). 



2.5. Processing 

After all maps have been convolved to SPIRE 350 /im 
resolution, we sample them with a hexagonal grid with 
approximately half-beam spacings (i.e. 12.5") . Uncer- 
tainties on the dust mass surface density, CO integrated 
intensity and HI column density are propagated through 
the necessary convolutions and samplings. Surface den- 
sities and other quantities have been deprojected using 
the orientation parameters listed in Table [I] 

3. SOLVING SIMULTANEOUSLY FOR a C o AND DGR 

In order to use the dust mass surface density to trace 
the total gas mass surface density, we assume 1) that 
dust and gas are well mixed (i.e. that Equation [3] holds), 
2) that the DGR is constant on ~kpc scales in our tar- 
get galaxies and 3) within a given ~kpc region, the 
DGR does not change between the atomic and molec- 
ular phases. Then, given multiple measurements of Ed, 
Shi and Ico that span a range of CO/HI values in a 
kpc-scale region of a galaxy, we can adjust aco till we 
find the value that returns the most uniform DGR for the 
region. This procedure makes no assumption about the 
value of the DGR, only that it is constant in that region. 
In addition, it makes no assumption about the value of 
aco o r the scales on which it varies. If aco varies on 
small scales, our measured aco value will represent the 
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Table 3 

Adopted Metallicities and Gradients 



PT05 


PT05 


KK04 


KK04 




Galaxy Central Metallicity 


Metallicity Gradient 


Central Metallicity 


Metallicity Gradient 


Source 


(12 + log(0/H)) 


(dex r^ 1 ) 


(12 + log(0/H)) 


(dex i" 1 ) 





NGC 0337 


8.18 


± 


0.07 






8.84 


± 


0.05 






M10 


Table 9 


NGC 0628 


8.43 


± 


0.02 


-0.25 


± 0.05 


9.19 


± 


0.02 


-0.54 


± 0.04 


M10 


Table 8 


NGC 0925 


8.32 


± 


0.01 


-0.21 


± 0.03 


8.91 


± 


0.01 


-0.43 


± 0.02 


M10 


Table 8 


NGC 2841 


8.72 


± 


0.12 


-0.54 


± 0.39 


9.34 


± 


0.07 


-0.36 


± 0.24 


M10 


Table 8 


NGC 2976 


8.36 


± 


0.06 






8.98 


± 


0.03 






M10 


Table 9 


NGC 3077 


8.30 


± 


0.20 






8.90 


± 


0.20 






Kll 




NGC 3184 


8.65 


± 


0.02 


-0.46 


± 0.06 


9.30 


± 


0.02 


-0.52 


± 0.05 


M10 


Table 8 


NGC 3198 


8.49 


± 


0.04 


-0.38 


± 0.11 


9.10 


± 


0.03 


-0.50 


± 0.08 


M10 


Table 8 


NGC 3351 


8.69 


± 


0.01 


-0.27 


± 0.04 


9.24 


± 


0.01 


-0.15 


± 0.03 


M10 


Table 8 


NGC 3521 


8.44 


± 


0.05 


-0.12 


± 0.25 


9.20 


± 


0.03 


-0.52 


± 0.15 


M10 


Table 8 


NGC S697 

1 N vJTVV OUZi I 


8.34 




0.24 






8.99 




0.10 






M10 




NGC 3938 


8.42 


± 


0.20 






9.06 


db 


0.20 






M10 


L-Z 


NGC 4236 


8.17 


± 


0.20 






8.74 


± 


0.20 






M10 


L-Z 


NGC 4254 


8.56 


± 


0.02 


-0.35 


± 0.08 


9.26 


± 


0.02 


-0.39 


± 0.06 


M10 


Table 8 


NGC 4321 


8.61 


± 


0.07 


-0.31 


± 0.17 


9.29 


db 


0.04 


-0.28 


± 0.11 


M10 


Table 8 


NGC 4536 


8.21 


± 


0.08 






9.00 


± 


0.04 






M10 


Table 9 


NGC 4569 


8.58 


± 


0.20 






9.26 


± 


0.20 






M10 


L-Z 


NGC 4625 


8.35 


± 


0.17 






9.05 


± 


0.07 






M10 


Table 9 


NGC 4631 


8.12 


± 


0.11 






8.75 


± 


0.09 






M10 


Table 9 


NGC 4725 


8.35 


± 


0.13 






9.10 


± 


0.08 






M10 


Table 9 


NGC 4736 


8.40 


± 


0.01 


-0.23 


± 0.12 


9.04 


± 


0.01 


-0.08 


± 0.10 


M10 


Table 8 


NGC 5055 


8.59 


± 


0.07 


-0.59 


± 0.27 


9.30 


± 


0.04 


-0.51 


± 0.17 


M10 


Table 8 


NGC 5457 


8.75 


± 


0.05 


-0.75 


± 0.06 


8.75 


± 


0.05 


-0.75 


± 0.07 


B07 




NGC 5713 


8.48 


± 


0.10 






9.08 


± 


0.03 






M10 


Table 9 


NGC 6946 


8.45 


± 


0.06 


-0.17 


± 0.15 


9.13 


± 


0.04 


-0.28 


± 0.10 


M10 


Table 8 


NGC 7331 


8.41 


± 


0.06 


-0.21 


± 0.31 


9.18 


± 


0.05 


-0.49 


± 0.25 


M10 


Table 8 



- jMoustakas et al.|j2010| M10), [KeTnicutt et al."]j2011| Kll), |Bresolin|(2007| B07). 



References. 



dust (or gas) mass surface density weighted average aco 
for the region. We discuss this process in detail in the 
following sections. 

A simultaneous solution for aco and DGR can only be 
performed if a range of CO/HI ratios are present in the 
measurements. The linear dependence of £h 2 on aco 
provides leverage to adjust aco m order to best describe 
all of the points with the same DGR over the range of 
CO/HI ratios. An "incorrect" aco value will cause a 
dependence of the measured DGR on CO/HI, increasing 
the spread in DGR values in the region. An illustration 
of this effect is shown in panel c) of Figure [I] Finding a 
solution or best-fit aco is equivalent to locating a mini- 
mum in the DGR scatter at a given value of aco- 

The basic procedure we use, minimizing the scatter 
in the PGR values in the region, was suggested in |Leroy| 
et al. (2011 ). There are, however, a variety of other tech- 
niques to solve for aco and DGR given multiple mea- 
surements, including directly fitting a plane to Ico? 
and Shi- It is not clear a priori which scatter mini- 
mization technique is optimal, so we performed a set of 
simulations, described in the Appendix, to test various 
techniques and optimize the procedure for our dataset 
and objectives. We describe the resulting scatter mini- 
mization procedure in more detail below. 

3.1. Defining the "Solution Pixel" 

To perform the solution, we require multiple measure- 
ments of dust and gas tracers in the region. We also aim, 
however, to select the smallest possible regions, in order 
to ensure an approximately constant DGR. We proceed 
by dividing each target galaxy into hexagonal regions en- 
compassing 37 of the half-beam spaced sampling points. 
We call these regions "solution pixels" (see panel a) of 



Figure [I] for an example). The 37-point solution pixels 
are a compromise between small region sizes and having 
a sufficient number of independent measurements to min- 
imize statistical noise. The solution pixels tile the galaxy 
with center-to-center spacing of 37.5". Thus, neighbor- 
ing solution pixels are not independent and share ~ 40% 
of their sampling points. The overlap between solution 
pixels is illustrated in Figure [T] Such a tiling is optimal 
because it fully samples the data and ensures that the 
final results are not dependent on the placement of the 
solution pixels. The solution pixels correspond to physi- 
cal scales ranging from ~0.6kpc to 3.9 kpc. 

3.2. Minimizing the DGR Scatter 

For each sampling point i in the solution pixel, the 
measurements of Ed,*, ^m,i an d Ico,i along with an 
assumed aco determine DGR^. We step through a grid 
of aco values to find the aco that results in the most 
uniform DGR for all DGR^ in the solution pixel. In all 
solutions presented here, we use an aco grid with 0.05 
dex spacing, spanning the range aco = 0.1 — 100 M 
pc -2 (K km s -1 ) -1 . 

We determine the most "uniform" DGR in a solution 
pixel by minimizing the scatter in the DGR^ values as a 
function of aco- The scatter is measured with a robust 
estimator of the standard deviation^] of the logarithm of 
the DGRi values — this technique appears to work best 
because outliers have little effect on the measured scat- 
ter because of both the logarithmic units and the outlier 
suppression. After measuring the scatter in the DGR at 
every aco value, we find the aco at which the scatter 

4 We use the IDL implementation of Tukey's biweight mean 
( |Press et al.|2002| biweightjnean.pro. 
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Figure 1. Illustration of the technique to determine DGR and aco f° r a single solution pixel in NGC 6946. Panel a) shows the CO 
map from HERACLES overlaid with the half-beam spaced sampling grid. The solution pixel in question is shown with a white hexagon 
and the 37 individual sampling points included in the solution are shown in red. A neighboring solution pixel is also highlighted with a 
thin white line to show how the pixels are arranged and that neighboring solution pixels share ~ 40% of their sampling points. Panel b) 
illustrates that the scatter in the DGR changes as a function of aco- Here we have plotted histograms of the measured DGR^ values in 
this solution pixel at three different values of aco (the optimal value in black and a factor of 5 above and below this value in green and 
blue, respectively). Panel c) illustrates that this scatter originates in variation of DGR as a function of CO/HI when aco is n °t at the 
optimal value. We show this by plotting the 37 DGR^ values as a function of CO/HI ratio for the same three aco values shown in panel 
b). The horizontal lines indicate the mean DGR for each set of points. The slope in DGR vs CO/HI space is minimized at the optimal 
aco- Finally, in panel d) we show the scatter in the DGR at each value of the full aco grid m black. The three highlighted aco values 
are marked with vertical lines. Panel d) highlights the fact that the DGR scatter is minimized at the best fit Log(aco)— 0.15 ±0.22 in this 
region. The minimization of the scatter in Log(DGR), as shown in panel d), is the technique we have determined to be the most effective 
for determining aco an d DGR, using tests that are described in detail in the Appendix. 



(Alog(DGR)) is minimized. This is taken to be our best 
fit aco value for the region. We consider a solution to 
be found when a minimum has been located in the DGR 
scatter within the range of our aco grid. This does not 
happen in every solution pixel — some do not have suffi- 
cient CO/HI contrast, others have too low S/N and for 
some the minimum is at the edge of the allowable range 
and the solution is not well-constrained. The failed solu- 
tions are not included in our further analysis. 

An illustration of the technique is shown in Figure [I] 
for a region in NGC 6946. In panel a) we show the 
HERACLES CO J=(2-l) map of the galaxy with our 



half-beam sampling grid, the hexagonal region shows the 
"solution pixel" in question, which includes 37 individual 
samples from the maps. Panels b) and c) show how vary- 
ing aco affects the mean DGR and scatter for the points 
in the region, illustrating how the scatter increases away 
from the best aco value. Panel d) shows the scatter as 
a function of aco f° r the whole aco grid. A clear mini- 
mum exists for this solution pixel at aco~ 1-4 M pc -2 
(Kkms- 1 )- 1 . 

3.3. Statistical Uncertainties on aco and DGR 
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We judge the uncertainties on the "best fit" aco and 
DGR in several ways. First, to take into account sta- 
tistical errors, we perform a Monte Carlo test on the 
solutions by adding random noise to our measured Ed, 
Ehi and Ico values according to each point's measure- 
ment errors. We repeat the solution with the randomly 
perturbed data values 100 times and find the standard 
deviation of the results. We also perform a "bootstrap- 
ping" trial, which tests the sensitivity of each solution 
to individual measurements. In each bootstrap iteration 
for a given solution pixel, we randomly select 37 sampling 
points, with replacement, and derive the solution. The 
bootstrap procedure is repeated 100 times for each solu- 
tion pixel and we measure the resulting standard devia- 
tion of the aco values. The standard deviations from the 
Monte Carlo and bootstrapping iterations are added in 
quadrature to produce the final quoted error for the aco 
values we determine. In addition, to check these uncer- 
tainties we also estimate scatter and bias in aco for the 
given technique from our simulated data trials described 
in the Appendix, based on the median CO S/N in the so- 
lution pixel and the measured minimum of Alog(DGR). 
The uncertainties from Monte Carlo plus bootstrapping 
are comparable to what we expect given the simulated 
data trials. 

3.4. Systematic Uncertainties on aco 

3.4.1. Uncertainties on aco from R21 Variations 

In the following, we report aco appropriate for the 
(1 — 0) line, since it is the canonical CO-to-H2 conver- 
sion factor that most observational and theoretical stud- 
ies utilize. To do so we have converted between (2 — 1) 
(which we have directly measured) and (1 — 0) using a 
fixed line ratio R21 = 0.7. This value of R21 is the average 
found for the entire HERACLES sample (E. Rosolowsky, 
private communication). Deviations from the R21 value 
we have assumed will result in systematic offsets in the 
(1 — 0) conversion factor, while the (2 — 1) conversion 
factor will be unaffected since it is what we have directly 
measured. To quantify the potential (1 — 0) aco offsets, 
we have investigated the variability of R21 in those galax- 
ies with publicly available CO J= (1 — 0) m aps from the 
Nobeyama survey of nearby spiral galaxies (Kuno et al. 
2007 ). Galaxies that have Nobeyama maps are marked 
with an asterisk in Table [IJ The details of this com- 
parison can be found in the Appendix. We find that 
deviations from R21 = 0.7 can cause small systematic 
shifts in aco appropriate for the (1 — 0) line, but the 
magnitude of the shifts are generally within the uncer- 
tainties on the aco solutions (i.e. typically less than 0.2 
dex). We note that variations of R21 within a pixel would 
introduce additional systematic uncertainties on aco- 

3.4.2. Variations of Ed Linearity Within Solution Pixels 

An assumption we make in using this technique is that 
the dust tracer we employ (here Ed) linearly tracks the 
true dust mass surface density. For the purposes of char- 
acterizing systematic uncertainties on aco, all that mat- 
ters is that within a given solution pixel Ed is a lin- 
ear tracer. This is because we calibrate DGR based on 
the Ed values themselves, so any multiplicative constant 
term cancels out in Equation [3] In theory, we could use 
a different linear dust mass surface density tracer in each 



pixel and still derive aco without issues. Systematic un- 
certainties in the absolute calibration of Ep s how up in 
the DGR and we discuss those in Section [531 

Non-linearity may be introduced in Ed within a so- 
lution pixel if dust emissivity or DGR vary with ISM 
phase. A variety of observations have suggested that dust 
emissivity increases in molecular gas relative to atomic 
gas (note, however, that most of the studies use CO to 
trace molecular gas and may interpret variations in aco 
as changes in emissivity). Recent work has suggested 
that the dust emissivity increases by a factor of ^ 2 be- 
tween the atomic and molecular ISM (iPlanck Collabo- 



ration et al. 2011b Martin et al. 20121). 11 the dust in 



molecular regions has a higher emissivity, Ed will over- 
estimate the amount of dust there, causing us to overes- 
timate the amount of gas. In that case, we would recover 
a higher aco than actually exists. As a first approxima- 
tion, our measured aco would be too high by the change 
in emissivity between atomic and molecular phases, a 
factor of ~ 2 based on the previously mentioned results. 
An increase in the DGR between atomic and molecular 
phases would also behave in this way, but we expect such 
ch anges to be less t han a factor of two (see Section 4.5 
of |Leroy et aLpOU) for more detailed arguments). 

Another potential source of non-linearity in Ed is ISM- 
phase dependent systematic biases in measuring Ed from 
SED modeling. Because warm dust will radiate more 
strongly per unit mass than cold dust at all wavelengths, 
the SED will not clearly reflect the presence of cold dust 
unless it dominates the mass. This means that the SED 
fitting technique is not sensitive to cold dust contained in 
GMC interiors (Ay > 1) at our spatial resolution. The 
fraction of the dust mass in these interiors, assuming a 
spherical cloud with uniform density and total Ay « 8 
mag, is ~ 40% (this estimate agrees well with recent ex- 
tinction mapping measurements of Milky Way GMCs: 



Kainulainen et al. 2011 Lombardi et al. 2011). If we 
underestimate the mass ol dust by missing cold dust in 
GMC interiors, we would underestimate the amount of 
molecular gas and adjust aco downwards. The magni- 
tude of this effect is at the factor of ~ 2 level and oppo- 
site in direction to what we expect for dust emissivity or 
DGR increases in molecular clouds. 

3.4.3. Opaque H I 

The HI maps we us e have not been cor rected for any 
optical depth effects (Walter et al. 2008). HI observa- 
tions of M 31 at high spatial and spectral resolution have 
found that the total gas m ass taking into acco unt opacity 
effects was ~ 30% higher (Bra un et al.||2009[ ). If optical 
depth effects are important, the optically-thin estimate 
for the atomic gas mass would be too low, resulting in 
our procedure determining a DGR that is too high (ex- 
cess dust compared to the amount of gas) . In the molec- 
ular regions, then, we will expect too much gas based on 
that same DGR, and consequently increase aco- There- 
fore, the effect of HI opacity is to overestimate aco- We 
expect the magnitude of such an effect to be a factor of 
~ 1.3 or less, which is well within the statistical uncer- 
tainties we achieve on the aco measurements. 

3.5. Systematic Uncertainties on DGR 
3.5.1. Absolute Calibration of Ed 
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As we have discussed above, as long as £d is a linear 
tracer of true dust mass surface density within a given 
solution pixel, its absolute calibration has no effect on 
the aco value we measure. The same is not true for 
the DGR value. Any uncertainties on the calibration of 
Ed will be directly reflected in the DGR meas urement. 
The Ed values we used are from fits of the |Draine fc] 



LB 02007) models to the IR SED using the Milky Way 
Tty = 3.1 grain model. The extent to which the appro- 
priate dust emissivity k v deviates from the value used 
by this model represents a systematic uncertainty on the 
DGR values we derive. Our knowledge of k v in different 
environments is limited, but there are constraints from 
observations of dust extinction curves and depletions in 
the Large and Small Magellani c Clouds (c.f. |Weingartner| 
fc Draine|200ip . |Draine et ah] Q2007| ) demonstrated that 
the Ud values decreased by a factor of ~ 1.2 when the 
LMC or SMC dust model was used instead of the MW 
model. Given that our sample is largely dominated by 
spiral galaxies and hence does not probe environments 
with metallicity comparable to the SMC (due to the 
faintness of CO in such regions and our S/N limitations), 
we expect that the systematic uncertaintie s on our DGR 
when comparing with other results from Dr aine fc Li| 
( |2007| ) model fits is small. However, other dust models 
may have different prescriptions for k v , so comparison 
of our PGR values to results from studies not using the 
Draine & Li (2007) models will have systematic offsets. 



4. RESULTS 
4.1. NGC 0628 Results Example 

We divided each of the 26 galaxies in our sample into 
solution pixels and performed the simultaneous solution 
for DGR and aco in each pixel. As an example, we 
present the results for NGC 0628 in the following section. 
The results for all solution pixels in all galaxies can be 
seen in the Appendix. 

Figure [2] shows, from left to right, the HI, CO and 
Ed maps used in our analysis. The circles overlaid on 
the map represent the centers of the solution pixels we 
have defined. Figure [3] shows the same circles represent- 
ing the solution pixel centers, now filled in with a color 
representing the best aco solution for that pixel on the 
left, a gray scale showing the uncertainty on that aco 
solution in the middle, and the DGR on the right. Fi- 
nally, in Figure II] we show these measured aco values as 
a function of galactocentric radius (^5). For comparison, 
Figure [3] also shows the local Milky Way aco= 4.4 M 
pc -2 (Iv km s -1 ) -1 with a solid horizontal line (dot ted 
lines show a factor of 2 above and below, see Section [5TT] 
for details on how the Milky Way value has been mea- 
sured). We note that the Milky Way may show a gra dient 
of aco with radius (also discussed in Section 5.1), but 
for purposes of comparison with the most widely-used 
conversion factor we use a constant aco on all plots. 

In general, NGC 0628 shows aco consistent with the 
Milky Way value within a factor of 2 at all galactocentric 
radii. Outside a radius of r2b ^ 0.6 we find few good 
solutions. There is a weak trend for lower aco a t smaller 
radii, with the central solution pixel having a conversion 
factor aco= 2.2 M pc -2 (K km s -1 ) -1 . Recent work 
by Blanc et al. (submitted) found consistent results for 
aco by inverting the star-formation rate surface density 



map using a fixed molecular gas depletion time. 

4.2. Completeness of Solutions 

The technique we have used to solve for DGR and aco 
simultaneously only works if there is sufficient S/N in the 
CO map and a range of CO/HI ratios in each solution 
pixel. These two constraints impose limits on where in 
the galaxies we can achieve solutions. In order to perform 
statistical tests on our sample of aco an d DGR values, 
we need to understand what biases these limits introduce 
into our results. For example, the failure of the technique 
in regions with low CO S/N generally limits our good 
solutions to the inner parts of galaxies, where metallicity 
tends to be higher. In order to judge the existence of 
trends in aco versus metallicity, we therefore need to 
understand where we have achieved good solutions. 

To investigate these effects, we have examined the frac- 
tion of solution pixels that have solutions and the uncer- 
tainty on the derived aco as a function of the range of 
CO/HI ratios and mean Ico in a pixel. In general, the HI 
in our target galaxies has a quite flat radial prof ile, while 
the CO dro ps off approximately exponentially (Schruba 
et al. 2011). Because of these trends, wherever there is 
sufficient signal in CO to achieve a good solution, the 
CO/HI range is adequate as well. Thus, we find that the 
mean Ico of a solution pixel is the best predictor for the 
existence and quality of a solution. Figure [5] shows the 
fraction of pixels with solutions and the average uncer- 
tainty on the solutions as a function of the mean Ico- 
We identify a cut-off at Ico> 1 K km s _1 above which 
we obtain solutions >90% of the time and those solu- 
tions have an average uncertainty of < 0.5 dex. For any 
statistical analysis that follows we use only pixels above 
this cut-off. 

4.3. Properties of Full-Sample aco Solutions 

In Tables [2] and [5] we list the average aco derived 
galaxy by galaxy and for the full sample. These averages 
only include solution pixels where Ico> 1 K km s _1 . The 
first two columns of Table [U list the mean and standard 
deviation of the individual solution pixel measurements. 
This approach treats each solution pixel equally, regard- 
less of how strongly it contributes to the total molecular 
gas mass. The third column instead lists the mean aco 
derived from the total £h 2 and the total Ico for the 
galaxy (above Ico=l K km s _1 ), which is equivalent to 
a mean where the pixels are weighted by their Ico- If 
a single aco value were to be applied to the data, this 
would be the optimal value to use. For most galaxies, the 
CO weighted mean is higher than the straight mean of 
the solution pixels. This is due to the fact that the molec- 
ular gas is not evenly distributed across the galaxy and 
the area- weighted average is different then the molecular 
gas weighted average. 

Figure [6] shows several histograms illustrating the dis- 
tribution of our measured aco values. The mean, de- 
rived with several different weighting schemes, is listed 
in Table [5] The top panels of Figure [6] show histograms 
of all solution pixels while the bottom panels show his- 
tograms of the galaxy averages from Table [4] On these 
histograms we highlight galaxies with high inclinations 
in green (i > 65°). It is clear that the high inclination 
galaxies tend to have higher aco on average than the 
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Figure 2. The HI, CO and Ed maps for NGC 0628, from left to right. The centers of the pixels in which we perform the simultaneous 
aco and DGR solutions are shown as circles overlaid on the images. The gray cross shows the central solution pixel for the galaxy. On 
the middle panel, the coverage of the HERACLES CO map is shown with a dotted line. Similar plots for all galaxies in the sample can be 
found in the Appendix. 
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Figure 3. Results of the simultaneous aco an d DGR solutions for NGC 0628. The centers of the solution pixels are represented with 
circles as shown on Figure [2] The left panel shows the resulting acO; the middle shows the uncertainty on that value and the right shows 
the DGR. Solution pixels where the technique failed do not appear. For NGC 0628 it is clear that where there are good solutions (as judged 
by the uncertainty on aco m the middle panel) most of the values are close to the MW Log(aico) — 0.64. In pixels with good solutions, 
the DGR varies smoothly across the galaxy, which shows that our assumption of a single DGR in each solution pixel is self-consistent. 



more face-on galaxies. In the highest inclination galax- 
ies the pixel will include contributions from gas at larger 
radii which tends to have lower DGR and less mole cular 
gas. This is equivalent to the challenge faced by |Leroy| 
et al. (2011) in the SMC, where some of the HI along 
the line of sight originates in an essentially dust-free en- 
velope, and in M 31, where regions on the minor axis 
of the galaxy have contributions from gas and dust at 
a variety of radii. In addition, optical depth effects for 
HI may be accentuated for highly inclined galaxies. All 
galaxies with i > 65° show average aco above the mean 
(except for NGC 0925, which ha s a somewhat uncertain 
inclination; de Blok et al. 2008). We thus eliminate all 
galaxies with i > 
ments. 



65 u , leaving 782 total aco measure- 
In Table [5] we list the average values of aco for the full 



galaxy sample. Excluding the high inclination galaxies, 
we find an average aco= 2.6 M pc -2 (K km s -1 ) -1 
for the individual solution pixel results. Weighted by 
Ico> the average value is slightly higher, aco= 2.9 M 
pc -2 (K km s -1 ) -1 . To avoid highly resolved galaxies 
like NGC 5457 and 6946 contributing more points to the 
average, we also calculate averages where each galaxy 
contributes uniformly. These are listed in the last two 
columns of Table [5] The average value for our sample, 
<^co= 3.1 M pc (K km s -1 ) -1 , is only slightly lower 
than what is found in the Milky Way disk. 

The standard deviation in aco is 0.38 dex, treating 
all lines of sight equally (with our Ico and inclination 
cut-offs). A key question we would like to answer is to 
what degree this scatter represents 1) true scatter within 
each galaxy, 2) galaxy-to-galaxy offsets or 3) variation 
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Figure 4. cnco solutions for NGC 0628 as a function of galactocentric radius (in units of r2s). The solid horizontal line shows the MW 
<^CO=4.4 Mq pc -2 (K km s -1 ) -1 (note that the MW could potentially have a gradient in aco with radius that we do not show here). 
The dotted lines show a factor of 2 above and below the MW value. The dashed horizontal line shows the average value for NGC 0628. 
The gray scale color of the points represents the uncertainty on aco as shown in Figure [3] — darker points have lower uncertainties. For 
comparison with the gray color table, two representative error bars for the aco solutions are shown in the top left corner of the plot. For 
NGC 0628, almost all of the high confidence aco solutions are within a factor of 2 of the MW value. 



Table 4 

Galaxy Average aco 



Galaxy Mean a C o Std. Dev. CO Weighted # Meas. 







(dex) 


Mean aco 




NGC 0337 


22.4 


16.2-31.0 (0.14) 


21.8 


2 


NGC 0628 


3.9 


2.1-7.4 (0.28) 


5.1 


67 


NGC 0925 


10.0 


6.8-14.7 (0.17) 


10.0 


1 


NGC 2841 


5.0 


2.5-10.3 (0.31) 


5.7 


19 


NGC 2976 


3.3 


0.9-12.5 (0.58) 


4.7 


18 


NGC 3077 


4.6 


2.3-9.4 (0.31) 


5.4 


7 


NGC 3184 


5.3 


2.9-9.5 (0.26) 


6.3 


43 


NGC 3198 


11.0 


6.4-18.9 (0.24) 


11.9 


11 


NGC 3351 


2.7 


1.0-6.9 (0.41) 


2.9 


28 


NGC 3521 


7.6 


4.9-11.8 (0.19) 


7.3 


42 


NGC 3627 


1.2 


0.4-3.3 (0.44) 


1.8 


43 


NGC 3938 


5.5 


4.1-7.5 (0.13) 


5.8 


19 


NGC 4236 









NGC 4254 


3.4 


2.1-5.7 (0.22) 


4.7 


46 


NGC 4321 


2.2 


1.1-4.6 (0.32) 


2.2 


57 


NGC 4536 


2.6 


1.0-6.7 (0.41) 


2.6 


13 


NGC 4569 


1.1 


0.3-4.1 (0.57) 


1.2 


14 


NGC 4625 









NGC 4631 


10.8 


5.6-19.5 (0.26) 


9.8 


40 


NGC 4725 


1.2 


0.4-3.2 (0.44) 


1.8 


7 


NGC 4736 


1.0 


0.5-2.0 (0.29) 


1.1 


33 


NGC 5055 


3.7 


1.9-7.4 (0.30) 


4.0 


86 


NGC 5457 


2.3 


1.1-4.8 (0.32) 


2.9 


142 


NGC 5713 


4.6 


1.7-12.6 (0.44) 


5.4 


13 


NGC 6946 


2.0 


0.9-4.4 (0.35) 


1.8 


158 


NGC 7331 


9.8 


6.2-15.3 (0.20) 


10.7 


32 


Note. — 


Averages include only solution pixels 


with Ico> 1 


K km s -1 . 



The number of measurements meeting this criterion are shown in the last 
column of the table. 
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Table 5 

Sample Average aco 



Sample 


Mean aco 


Std. Dev. 


CO Weighted 


Gal. Weighted 


Gal & CO Weighted 




all L.O.S. 


(dex) 


Mean aco 


Mean aco 


Mean aco 


Incl < 65° 


2.6 


1.0-6.6 (0.41) 


2.9 


3.1 


3.5 


Incl > 65° 


7.2 


3.0-17.6 (0.39) 


8.2 


6.5 


6.7 



Note. — Averages include only solution pixels with Ico> 1 K km s 
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Figure 5. The completeness and quality of our aco solutions as a 
function of the mean Ico m a solution pixel. The black points show 
the fraction of solution pixels where a solution was found. The gray 
points show the mean uncertainty on aco with error bars showing 
the standard deviation. For Ico> 1 K km s _1 , > 90% of solution 
pixels have solutions and their typical uncertainty is < 0.5 dex. 
We use this 90% completeness cut to select a well-defined sample 
of aco measurements for statistical analysis in the following work. 

of aco as a function of local parameters. If the scatter 
comes from variations with local environmental parame- 
ters, we may be able to generate a prescription for aco 
as a function of other observables. In the following sec- 
tions we explore the variations of aco within and among 
our galaxies to understand if and why it varies. 

4.4. Radial Variations in aco 

In Figure [7] we present a summary of the aco val- 
ues we find as a function of galactocentric radius. Plots 
of aco from each individual galaxy as a function of T25 
are presented in the Appendix. Each individual solution 
that makes our Ico an d inclination cut is shown in Fig- 
ure [7] with gray circles. We also display the mean and 
standard deviation of the aco values in 0.1 T25 bins for 
each galaxy (this mean treats all solution pixels equally) . 
The top panel of Figure [7] shows the measured aco and 
the bottom panel shows those same values normalized by 
each galaxy's average aco from Table |4j 

The average radial profile of our galaxies is mostly flat 
as a function of r25, but almost all galaxies show a de- 
crease in aco m the inner ~ 0.2 r 2 5 compared to their 
average value. The mean central decrease is ~ 0.3 — 0.4 
dex, but it can be as much as 0.8 dex. For several galaxies 
this decrease leads to a central aco value that is an order- 
of-magnitude lower than the MW aco (e.g. NGC 4736, 
5457 and 6946). The central depression persists even 
after correcting for differences in R21 (see the following 
Section for further discussion). 

Normalizing each galaxy by its average aco and mak- 
ing each galaxy contribute equally to the average leads 
to much smaller scatter at T25 > 0.3 (this can be seen 
in the bottom panel of Figure [7] where the radial pro- 
files show much less scatter than the individual measure- 
ments, which is not the case in the top, unnormalized 



panel). Inside that radius, scatter in the normalized aco 
profiles is reduced, but not by as much outside that ra- 
dius, indicating that individual galaxies show different 
central profiles of aco- The overall conclusions from ex- 
amination of Figure u\ are: 1) outside Y25 ^ 0.3 most of 
the scatter in our aco measurements can be explained 
by galaxy to galaxy differences, 2) inside that galacto- 
centric radius, galaxy-specific trends in aco dominate 
the scatter and 3) in general the inner region of galaxies 
shows aco lower than the rest of the disk. 

4.5. aco in Galaxy Centers 

In Table [6] we list the central aco values and compare 
them to the MW value and the mean for the galaxy. Al- 
most all of the central solution pixels have aco lower 
than the galaxy mean. The galaxies NGC 3351, 3627, 
4321, 4736, 5457 and 6946 show aco values a factor of 
5 or more lower than the MW value with 3<r confidence. 
For 4736, 5457 and 6946, in particular, the central solu- 
tion pixel has aco an order-of-magnitude below the MW 
value. Only one galaxy, NGC 0337, shows a central aco 
higher than the MW value at 3a. In Table [6| we also list 
the CO (2 — !)/(! — 0) ratio measured as described in 
Section [3.4. 1| While the central solution pixels do tend 
to have higher R21 than the value we have adopted, the 
difference is less than a factor of 2 in all cases, and does 
not alter our main conclusions about the low (1 — 0) con- 
version factor in these regions. In particular, for the 10 
galaxies with measured R21 and i < 65°, the central aco 
corrected for R21 is still on average 0.3 dex lower (i.e. a 
factor of 2) than the galaxy average. 

The central solution pixels can be outliers from the 
rest of the pixels in our sample in the sense of having 
high CO/HI ratios. Since there are relatively few so- 
lution pixels with these conditions, they may not have 
been well represented in the simulated data we used to 
test the accuracy of the solution technique. High CO/HI 
ratios could potentially bias the results: if the S/N of the 
HI maps is almost always higher than the S/N of CO, 
the DGR scatter could be reduced in these conditions 
purely by decreasing the importance of CO in assessing 
the gas mass surface density. Since this bias would move 
the aco results towards lower values in the centers, we 
performed a test to judge whether our low central aco 
measurements could be due to this effect. The details of 
the test are described in the Appendix. In brief, we gen- 
erated simulated datasets with known aco values where 
the Nhi, Ico an d S/N is matched to the observations 
in each central solution pixel and performed Monte Carlo 
trials to see how well the known input aco was recov- 
ered. In all cases, we found no evidence for a bias in the 
central solution pixel. 

Taking into account the offset from R21 variations, our 
results show that the central aco m NGC 3351, 3627, 
4321, 4736, 5457 and 6946 is lower than the MW value 
by a factor of 4 — 10. Many of these galaxies show aco 
closer to the Milky Way value at larger radii, in line with 
the general trend seen in Figure [7) For NGC 3627, the 
mean aco for the galaxy is low as well, suggesting that 
the central region is not distinct from the rest of the disk 
(this galaxy is an interacting member of the Leo Triplet, 
so it is unique in our sample). Conversely, examination 
of the individual galaxy aco measurements as a function 
of radius (shown in the Appendix) for NGC 3351, 4321, 
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Figure 6. Histograms showing the distribution of aco results with various weighting schemes. The averages are reported in Table [5] 
The panels show histograms of galaxies with i < 65° in black and galaxies with i > 65° in green. In each panel we highlight the mean 
value of aco f° r galaxies with i < 65° with a dotted vertical red line and the Milky Way aco m g ra y- These histograms include solution 
pixels above the Ico > 1 K km s — 1 cutoff. The top left panel shows the histogram of the solution pixels. The top right panel shows those 
same values weighted by their Ico — this makes pixels contribute to this histogram in proportion to their molecular gas mass. The bottom 
panels show histograms of the average aco from each galaxy as listed in Table PI Galaxies with high inclinations tend to show a higher 
than average aco, because the solution pixels include gas and dust at a range of radii. On average, we find aco slightly lower than the 
MW value in our sample. 



Table 6 

Central aco Measurements 



Galaxy 


Logio(a C 0,Cen) 

±a 


Amw = 
Logio(a C o/«CO,Mw) 




^mean — 

Logio(a C o/(«CO>) 


^mean/ & 


T21 I 


^ogi (r 2 i/0.7) 


NGC 0337 


+1.25 


± 0.17 


+0.61 


3.6 


-0.10 


0.6 






NGC 0628 


+0.35 


± 0.24 


-0.29 


1.2 


-0.24 


1.0 






NGC 2976 


+0.00 


± 0.49 


-0.64 


1.3 


-0.51 


1.0 






NGC 3077 


+0.60 


± 0.17 


-0.04 


0.2 


-0.06 


0.4 






NGC 3184 


+0.25 


± 0.16 


-0.39 


2.5 


-0.47 


3.0 


0.77 ± 0.05 


+0.04 


NGC 3351 


-0.15 


± 0.14 


-0.79 


5.6 


-0.58 


4.1 


1.10 ± 0.24 


+0.19 


NGC 3627 


-0.25 


± 0.14 


-0.89 


6.3 


-0.34 


2.4 


0.54 ± 0.05 


-0.12 


NGC 3938 


+0.70 


± 0.19 


+0.06 


0.3 


-0.04 


0.2 






NGC 4254 


+0.95 


± 0.82 


+0.31 


0.4 


+0.41 


0.5 


1.03 ± 0.06 


+0.17 


NGC 4321 


-0.20 


± 0.17 


-0.84 


4.8 


-0.55 


3.1 


1.25 ± 0.16 


+0.25 


NGC 4536 


+0.35 


± 0.13 


-0.29 


2.2 


-0.07 


0.5 


1.26 ± 0.28 


+0.26 


NGC 4625 


+1.05 


± 0.41 


+0.41 


1.0 










NGC 4725 


-0.15 


± 0.71 


-0.79 


1.1 


-0.22 


0.3 






NGC 4736 


-0.55 


± 0.17 


-1.19 


6.9 


-0.56 


3.2 


1.35 ± 0.09 


+0.29 


NGC 5055 


+0.00 


± 0.25 


-0.64 


2.6 


-0.57 


2.3 


1.10 ± 0.08 


+0.20 


NGC 5457 


-0.45 


± 0.20 


-1.09 


5.5 


-0.80 


4.1 


0.90 ± 0.08 


+0.11 


NGC 5713 


+0.50 


± 0.25 


-0.14 


0.6 


-0.16 


0.6 






NGC 6946 


-0.40 


± 0.31 


-1.04 


3.4 


-0.70 


2.3 


1.07 ± 0.14 


+0.19 



Note. — Galaxies with i > 65° have been omitted from the Table. 
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NGC4321 I 
NGC4536 I 
NGC4725 - 
6 - 



NGC5713 I 
NGC6946 : 
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Figure 7. Solutions for aco plotted versus their galactocentric radius (r2s). The top panel shows the aco solutions and the bottom panel 
shows the same values normalized by each galaxy's average aco (listed in Table (41). All individual solutions are shown as gray points. The 
mean and standard deviation of all of the solutions in 0.1 T25 bins are shown witholack symbols (this mean treats all lines of sight equally). 
The solid black horizontal line in the top panel shows aco f° r the Milky Way. The dotted black horizontal line shows the average, with no 
weighting, of all solution pixels: aico=2.6 Mq pc -2 (K km s -1 ) -1 . The radial profile for each galaxy is shown with a solid colored line. 
The average aco radial profile of our galaxies is mostly flat as a function of 125 with a decrease in aco towards the center of the galaxy. 
On average, the central aco value is ~0.3 dex lower than the rest of the galaxy. 
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4736, 5457 and 6946 demonstrates that these galaxies 
show an unresolved region of lower than average aco 
in their centers (note that our solution pixel grid over- 
samples the data, so an unresolved central depression 
affects the 125 = point and the six adjacent solution 
pixels). The closest galaxy that shows a central depres- 
sion is NGC 4736 at D=4.66 Mpc. At this distance, our 
solution pixel covers a region of radius 0.8 kpc. Even 
for this nearby example, we do not resolve the central 
depression. We tested whether the depressions were re- 
solved using a grid of independent solution pixels (i.e. 
not overlapping) and found that the depression only af- 
fected the central pixel, consistent with it being unre- 
solved. We note that our ability to detect any central 
depression may be a function of the galaxy's distance 
due to the increased size of the central solution pixel. 
Even at D= 14.3 Mpc, however, we detect a clear cen- 
tral depression in NGC 4321. 

The type of nuclear activity in each galaxy is not a 
good predictor for whether or not it displays a central 
aco depression. Of the galaxies that show the clearest 
central depression NGC 3351, 5457 and 6946 are clas- 
sified as star-formation or HII region dominated, while 
NGC 4321 and 4736 have signatures of active galactic 
nucleus (AGN) or low-ionization nuclear emission-line re- 
gion (LINER) act ivity (for details on the nuclear classifi- 



Table 7 

aco Correlation Properties 



cations, see |Kennicutt et aLpQll[ ). Several galaxies with 
flat radial aco profiles do show evidence for AGN activ- 
ity: NGC 3627, 4254, and 4725, for example. Enhanced 
AGN activity could poten tially effect mole cular gas prop- 



erties in the nuclei (e.g. Krips et aL||2008 ), but the rela- 
tively weak AGN present in the KIN GF1SH galaxies may 
not dominate on the kiloparsec scales we study here. 

To summarize, we find that several galaxies in our 
sample show central aco values that are substantially 
lower than the Milky Way value and also well below the 
galaxy average. This appears to take the form of de- 
pression in the central region that is unresolved by our 
solution pixel grid. We discuss these central regions and 
the physic al co nditions that may lead to low aco further 
in Section [631 



4.6. Correlations of aco with Environmental 
Parameters 

Variations in aco could potentially be related to 
variations in the local environmental conditions includ- 
ing metallicity, ISM pressure, interstellar radiation field 
strength, gas temperature, dust properties or other vari- 
ables. Correlations between aco and quantities that 
trace these environmental conditions may allow us to 
identify the drivers of aco variations. They may also 
provide tools to predict the appropriate aco in a given 
environment. In the following, we examine the correla- 
tions of our measured aco with several observables that 
trace ISM conditions. For each tracer, we use the average 
value in each solution pixel. The physical interpretation 
of t hese correlations or lack thereof is discussed in Sec- 
tion ED 

In Figure [8] we plot our measured aco values as a func- 
tion of the average radiation field intensity (U), the PAH 
fraction (<?pah), metallicity (12+Log(0/H)), stellar mass 
surface density (£*), star- formation rate surface density 
(^sfr) and the dust mass surface density (Ed). Each 
panel shows the individual measurements with gray sym- 



Log(aco) vs. 
Variable 


Correlation 21 
Coeff. 


a from 
Null 


Linear Fit b 
Offset Slope 




+0.19 


5.2 


0.29 


+0.29 


Log(Z7) 


-0.19 


5.2 


0.50 


-0.33 


qpAH (%) 


+0.13 


3.7 


0.11 


+0.08 


12+Log(0/H) 


-0.05 


1.4 


2.39 


-0.24 


Log(£*) 


-0.25 


7.1 


1.05 


-0.33 


Log(E D ) 


-0.15 


4.1 


0.26 


-0.22 


Log(EsFR) 


-0.11 


3.2 


0.06 


-0.14 



Note. — These correlations are shown in Figure [8| 
a Spearman rank correlation coefficient. 

k Correlation with metallicity using a homogeneous sample of 
metallicity measurements is discussed in Section |4.6.1| 

bols and the average for the galaxies with red symbols. 
The Milky Way aco is highlighted with a horizontal gray 
line and the mean aco treating all solution pixels equally 
(<^co= 2.6 M pc -2 (K km s -1 ) -1 ) is shown with a 
dashed horizontal line. Overlaid on each panel is the 
binned mean and standard deviation and a linear fit to 
the gray points. In Table [7] we list the linear fit results 
as well as the rank correlation coefficient for each panel 
and its significance (in standard deviations away from 
the null hypothesis). The rank correlation coefficients 
suggest there are significant correlations between aco 
and all of the variables except metallicity. Due to poten- 
tial inconsistencies between various metallicity measure- 
ments and incomplete knowledge of metallicity gradients 
in some galaxies, we examine the trends with metallicity 



separately in Section 4.6.1 



Several other radiation field properties are measured 
from the dust SED modeling. These include U m i n , the 
minimum radiation field heating the dust, and fpDR the 
fraction of the dust luminosity that arises in "PDR-like" 
regions where U>100 Ummp- We have investigated the 
dependence of aco on these quantities and find weak 
trends with low significance. 

The existence of a correlation between aco and en- 
vironmental parameters does not directly identify the 
cause of the variations in aco, particularly since all of 
the parameters change radially to first order. In Figure [9] 
we illustrate the radial variations of the same variables 
by plotting them normalized by their average value as 
a function of r25- All of the variables show significant 
correlations with radius. Table |8] lists the correlation co- 
efficients and linear fits to the normalized radial profiles 
of the parameters. In all cases, the correlation of the vari- 
ables with r25 is more significant than the correlation of 
aco with those variables. 

If all of the variations of these parameters were pri- 
marily related to radius, we would not expect to find 
stronger correlations between aco and any parameter 
than between aco and radius. Our results suggest, how- 
ever, that the correlation between aco and is stronger 
(r s = —0.25) and more significant (7.1<j) than the corre- 
lation with radius Y25 (r s = +0.19, 5.2a). One possible 
explanation is that our measurement may trace the 
stellar profile better than T25 itself. Other potential ex- 
planations for the corre lation of aco with will be 
discussed in Section l6Jl 
Figure [7] illustrates that the normalized radial profile of 
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Figure 8. aco as a function of environmental paramet ers. The panels show aco plotted versus average radiation field intensity ([/; 
top-left) and PAH fraction (qpAH5 top-right) from the |Aniano et aL] (|2012) fits, metallicity (12+log(0/H); middle-left), stellar mass surface 
density (E*; middle-right), star- formation rate surface density (2Jsfr; lower-left) and dust mass surface density (Ed; lower-right). The 
solid gray line shows the Milky Way aco an d the dotted black line shows the sample average aco= 2.6 Mq pc -2 (K km s -1 ) -1 , treating 
all solution pixels equally. Individual aco solutions above our Ico an d inclination cuts are shown with gray symbols. The mean and 
standard deviation of those values in bins are shown with black circles and error bars. A linear fit to the measurements is shown with a 
black line. The mean values for all galaxies are shown with red circles. All variables except for metallicity show a statistically significant 
correlation with aco- None of these correlations allow aco to be predicted with significantly better precision than the scaling with T25. 
However, the extreme values of SgFR (> 0.5 Mq yr _1 kpc -2 ) and and £* (> 1000 Mq pc -2 ) are always associated with low aco an d 
are located in galaxy centers. 



aco shows a factor of ~ 2 standard deviation in a typical 
radial bin. This uncertainty can be explained primar- 
ily by the uncertainty on our aco solutions themselves, 
which is typically close to a factor of ~ 2. It is therefore 
unlikely that we could predict aco as a function of other 
variables to better precision than a factor of ~ 2 unless 
the galaxy average aco is highly correlated with that 
variable. In general, the galaxy averages shown in Fig- 
ure [8] do not appear to be more tightly correlated with 
the environmental parameters. Figure [8] thus illustrates 
that aside from outliers at the extremes of these plots, 
the minimum standard deviation of aco m these bins is 



a factor of 2 or more and correlations with environmen- 
tal parameters do not allow us to predict the behavior 
of aco within the galaxies better than our normalized 
radial profile. 

As previously discussed, the average profile of aco ver- 
sus radius is mostly flat with a central depression in some 
galaxies. Many of the environmental parameters we plot 
in Figures [8] and [9] show radial trends extending over the 
entire range we cover, although the normalized trends 
generally span less than an order of magnitude. This 
range does not provide the leverage to separate the dom- 
inantly radial correlations of multiple variables. Future 
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Table 8 

Correlation Properties vs. T25 



Figure 9. Normalized environmental parameters as a function of galactocentric radius r25- In each panel we show the average value of 
each quantity in the solution pixel, normalized by the mean value for the galaxy, with gray circles. The mean and standard deviation of the 
binned, normalized values are shown with black circles and error bars. A linear fit is shown with a solid black line. All parameters (radiation 
field strength, PAH fraction, metallicity, stellar mass surface density, star-formation rate surface density, dust mass surface density and aco) 
show significant radial correlations. Since we have assumed radial gradients in the metallicity for many galaxies, normalized metallicity 
directly reflects the slopes of these gradients. In general, these radial trends are much stronger than the correlations between the variables 
and aco, presented in Figure [8] 

studies attempting to associate changes in aco with en- 
vironmental parameters will need either higher precision 
measurements of aco or to span a greater range of envi- 
ronments. 

One possible way to overcome the limitations of sep- 
arating various radial trends is to normalize all of the 
variables by their mean in radial bins and then search for 
residual trends between them. We have investigated such 
non-radial variations in 0.1 T25 bins and in general, find 
only very weak trends. NGC 4254 is one of few galax- 
ies that shows a marginally significant correlation — the 
normalized aco measurements correlate with normalized 
Note. These correlations are shown in Figure [9] U and at the level of 4-5a from the Spearman rank 

Spearman rank correlation coefficient. correlation coefficient. In both cases, the correlation is 
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positive, such that aco increases in region with higher 
U and E* compared to the average in that radial bin. 
NGC 6946 also shows a correlation at 5<r between nor- 
malized ftco and fpDR- In this case, the correlation is 
negative, meaning aco decreases in regions with high 
fpDR- 

NGC 6946 is unique in that the non-radial structure 
is clearly visible in the maps of aco shown in the Ap- 
pendix. The low aco values in that galaxy appear to 
track the spiral arm structure seen in the Ico and Ed 
maps. To make this observation quantitative, we have 
investigated the correlation between the radially normal- 
ized aco, Ico, and Ehi- These correlations are weak 
(all are < 5<r signficance in the Spearman rank correla- 
tion) but more widespread. NGC 3627, 4254, 4321, 5457 
and 6946 all show some degree of correlation at > 3a 
with either radially normalized Ico, ^hi or Ed- For all 
galaxies aside from NGC 6946, the correlations are pos- 
itive, in the sense that the conversion factor is higher 
where there is more HI and dust or more CO emission. 
NGC 6946 shows negative correlations with all of these 
parameters. 

4.6.1. aco versus Metallicity 

Metallicity has been suggested by several theoreti- 
cal and observational studies to be an important driver 
for aco variations. Unfortunately, metallicity measure- 
ments in nearby galaxies are often very uncertain and 
subject to systematic errors from different calibrations 
and techniques. In Figure [8] we show our aco measure- 
ments as a function of metallicity from the PT05 calibra- 
tion (middle- left panel). Our measurements span ~ 0.5 
dex in metallicity with no statistically significant trend 
in aco- One possible reason for the lack of a clear trend 
in this figure is that we have combined metallicity mea- 
surements obtained with different techniques. In addi- 
tion, many of the galaxies in our sample lack constraints 
on possible gradients. 

To explore any metallicity trends that may be washed 
out due to systematic effects when combining metallici- 
ties from different sources, we isolate a sample of galaxies 
with uniform ly determined metallicities from Moustakas 



et al. (2010) and plot those separately in Figure 10 We 
have eliminated galaxies whose metallicities have been 
determined from integrated spectrophotometry (those 
listed as "M10 Table 9" in Table m because the drift 
scan observations were not always along angles that al- 
low a robust gradient measurement. We also isolate our 
highest confidence aco measurements by showing only 
those with uncertainties less than 0.3 dex (a factor of 
- 2). We show metallicities in both the PT05 and KK04 
calibrations for comparison. 

In general, there is a slight negative correlation of aco 
with metallicity for this subsample. In the PT05 calibra- 
tion, we find a rank correlation coefficient of r = —0.2, 
which is 3.8<j from the null result given the number 
of measurements. Using KK04 metallicities, we find 
r = —0.1 at 2.0<j from null. If we remove the S/N cut on 
the aco measurements, the correlations in both PT05 
and KK04 essentially disappear, yielding results that are 
less than la from the null result. 

The weakness of the correlation between aco and 
metallicity suggest that in the regions of the galaxies we 
are studying, metallicity may not be the primary driver 



of aco variations. Along these lines, it is interesting to 
note the contrast between a galaxy like NGC 6946, which 
has a relatively shallow metallicity gradient, and galaxies 
like NGC 0628 and 3184 which are thought to have much 
steeper gradients. NGC 6946 has aco spanning a range 
of an order of magnitude, while NGC 0628 and 3184 have 
much smaller ranges — the opposite of what we would ex- 
pect if their metallicity gradient was the dominant factor 
controlling aco- 

4.7. Properties of Full-Sample DGR Solutions 

In the Appendix we show the measured DGR for all 
galaxies and solution pixels in figures similar to Figure [3j 
The DGR we report is the mean value for all of the indi- 
vidual sampling points in a solution pixel. Generally, the 
scatter in the DGR^ values within a solution pixel is small 
(< 0.1 dex). This means that the uncertainty on the 
aco value generally dominates the uncertainty on DGR 
as well. Therefore, it is important to note that the errors 
in the DGR and aco values are highly correlated. We 
assign a representative uncertainty on the DGR using the 
±lcr bounds on aco- In general, when good solutions are 
obtained, the DGR measurements vary smoothly across 
the galaxy. This supports our key assumption that DGR 
varies on scales larger than our solution pixels. The 
smoothness of the DGR maps can be seen by inspecting 
the figures in the Appendix and the results for NGC 0628 
shown in Figure [3J 

In Figure 11 we show a histogram of all of the measured 
DGR valueslor our galaxies that comply with our cuts 
on inclination and Ico- The average of our measure- 
ments treating all lines of sight equally is Log(DGR) = 

— 1.86 with a standard deviation of 0.22 dex. This is 
slightly higher than the value typically adopted for the 
solar neighborhood of Log(DGRMw)= —2.0. Forcing 
all galaxies to contribute equally to the average despite 
differing numbers of solution pixels gives Log(DGR) = 

— 1.96. The lower scatter in the DGR measurements for 
all solution pixels compared to the aco measurements, 
despite the uncertainty on aco dominating the uncer- 
tainty on DGR, is due to the contribution of HI to the 
total gas mass surface density. Because HI makes up 
some fraction of the gas mass the DGR has a smaller 
possible range over which to vary compared to aco- 

4.8. Correlations of DGR with Environmental 

Parameters and Radius 

DGR may also vary as a function of environmental 
parameters. To explore any such trends we show our 
measured DGR values as a function of the same tracers 



we have previously studied in Figure [12] We list the cor- 
relation coefficients and linear fit parameters in Table [9j 

DGR shows significant correlations with galactocentric 
radius, metallicity and E*. Even with heterogeneously 
determined metallicities, it shows a much clearer trend 
with 12+log(0/H) than aco does over the same range. 
A correlation of DGR with metallicity is expected — given 
the high depletions of elements such as Mg, Si, Ca and 
Ti in the local area of the Milky Way, to first order the 
mass of dust should be proportional to the amount of 
heavy elements. 

In Figure [l3j we show DGR as a function of metallicity 
for the same sample of galaxies from Figure [To] In addi- 
tion to the best linear fit, shown with a dotted black line, 
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Figure 10. Measured aco values as a function of metallicity for galaxies with HII region abundance gradients from Moustakas et al.| 
(2010) in the PT05 (left) and KK04 (right) calibrations. Only measurements with uncertainties less than < 0.3 dex are shown. The mean 
and standard deviation in 0.1 dex bins of metallicity are shown with yellow circles and black error bars. A linear fit to aco(Z) is overlaid 
with a solid black line. The correlation of aco with metallicity is weak for both calibrations. In particular it is clear that galaxies with 
shallow metallicity gradients like NGC 6946, can still span a wide range of aco values, whereas those with large gradients, like NGC 3184, 
may have essentially constant aco- 



All Solution Pixels 



Galaxy Averages 



200 



150 - 



100 - 




-2,5 -2,0 -L5 
Log(DGR) 



L0 -3,0 -2,5 -2,0 -1,5 
Log(DGR) 



-1,0 



Table 9 

DGR Correlation Properties 



Figure 11. Histograms of the measured DGR values for all solution pixels (left) and galaxies (right) in our sample. The dotted line shows 
the mean value for the sample and the red line shows the local MW DGR of 0.01 for comparison. The average DGR for our sample is 
slightly higher than in the MW. The standard deviation of the DGR measurements is 0.22 dex, which is smaller than the scatter in the 
aco values for these same pixels. 

that some dust is associated with HI. 

The DGR values show a stronger correlation with 
metallicity than the aco values. In the PT05 calibra- 
tion we find a rank correlation coefficient of 0.35 which 
is significant at 6.4<j. For KK04, the correlation coef- 
ficient is 0.39 at 6.9a. In a given metallicity bin, the 
standard deviation of the DGR values is ~ 0.15 dex for 
PT05 and - 0.18 for KK04. This suggests that, on av- 
erage, we should be able to predict DGR to better than 
a factor of 2 given the metallicity in one of these cali- 
brations. In addition, we do not see evidence for major 
galaxy-to-galaxy offsets in the DGR versus metallicity 
plots, which distinguishes the DGR and aco behavior, 
and most galaxies appear to have similar slopes in the 
plots. 

Neither the KK04 nor PT05 calibration values fall di- 
rectly along the line of scaled MW DGR. On the PT05 
scale, almost all of our sample has higher DGR for a given 
metallicity than the Milky Way scaling, while the oppo- 
site is true on the KK04 scale. For purposes of scaling 
the local Milky Way metallicity, we have used strong-line 
abundance measurements calculated from the integrated 
spectrum of the Orion Nebula (Croxall et al, in prep). 
The best linear fit to the data is shown with a dot- 
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Note. — These correlations are shown in Figure |8| 
a Spearman rank correlation coefficient. 

k Correlation with metallicity using a homogeneous sample of 
metallicity measurements is discussed in Section |4.8| 



we also plot a prediction for scaling the local MW DGR 
linearly with metallicity (solid black line) and a factor of 
2 higher and lower DGR with dashed lines. The scaling 
of the plot is such that it covers three orders of magni- 
tude, the same range covered in Figure [lOj It is clear 
that the DGR values have much less scatter than the 
aco values over the same range of metallicity — a prod- 
uct of the limited allowable range of DGR set by the fact 
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Figure 12. DGR as a function of environmental paramet ers. The panels show DGR plotted versus average radiation field intensity (U; 
top-left) and PAH fraction (qpAH 5 top-right) from the |Aniano et al.| (|2012) fits, metallicity (12+log(0/H); middle-left), stellar mass surface 
density (£*; middle-right), star- format ion rate surface density (^sfr; lower-left) and galactocentric radius (r2s; lower-right). The solid 
gray line shows the Milky Way DGR and the dotted black line shows the sample average. Each individual solution above our Ico an d 
inclination cuts are shown with gray symbols. The mean and standard deviation of those values in bins are shown with black circles and 
error bars. A linear fit to the measurements is shown with a black line. The mean values for all galaxies are shown with red circles. DGR 
shows statistically significant correlations with r25, 12+log(0/H) and £*. 



ted line in both the PT05 and KK04 panels. This fit is 
parametrized by the following equation where uncertain- 
ties have been determined using bootstrapping: 



Log(DGR) 



6(12 + Log(0/H) 



(4) 



For the PT05 measurements we find the following con- 
stants: 

a = -1.86 ±0.01 
b = 0.85 ±0.11 
c = 8.39 

For the KK04 measurements the constants are: 
a = -1.86 ±0.01 



b = 0.87 ±0.11 
c = 9.05 

In these equations, the constant c is the mean metallicity 
of the sample of points. We use this parametrization in 
order to avoid covariance in the offset (a) and slope (b). 
In both cases the slope is below unity. It can be seen 
in Figure [13] that the flatter than unity slope is mainly 
due to NGC 3184 (blue stars). Leaving this galaxy out 
of the fits produces slopes of 1.13±0.10 and 1.01±0.11 
for PT05 and KK04, respectively. It also increases the 
significance of both correlations by l-2a. 

5. COMPARISON TO LITERATURE 
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Figure 13. Measured DGR values as a function of metallicity for galaxies with HII region abundance gradients from JMoustakas et al.| 
(2010l in the PT05 (left) and KK04 (right) calibrations. Only solutions where the uncertainty in aco is l ess than 0.6 dex are shown, 
although the significance of the correlation between DGR and Z is relatively unchanged if all points are included. A linear fit to DGR(Z) is 
overlaid with a dotted black line. There is a clear correlation of DGR with Z. The best fit slope is slightly below linear, but if NGC 3184 is 
removed the correlation is consistent with linear. The solid black line shows a linear scaling of the MW DGR and metallicity, with dashed 
lines above and below showing a factor of 2 difference from MW. On the PT05 scale, our sample appears to have higher DGR for a given 
metallicity compared to the MW, while the opposite is true on the KK04 scale. 



5.1. Measurements of aco in the Milky Way and 
Nearby Galaxies 

The local region of the Milky Way is seen to have a con- 
version factor close to aco = 4.4 M pc -2 (K km s -1 ) -1 . 
This standard value has been recovered to within a fac- 
tor of 2 by various techniques: 7 -ray measurem ents ([DTI 
gel et al.|1996 [ Abdo et al.|2010|); v irial masses (|Solomon| 



et al.||1987D7~and dust ([Dame et al.||2 001 
2008; Planck Collaboration et al.||201 



TaT 

rTnc 



Pineda et al.j 
A number ol 

studies ol the Galactic Center, however, ha ve found very 



different answers. Using 7-ray observations, Strong et al 
(2004) found evidence that aco m the Galactic Center 



was factors of 5—10 lower than in the disk. Similar results 
using dust as a tracer for total gas mass have been found 
by |Sodroski et al.| ( |1995| ) who suggest that aco is lower 
by factors of 3-1 Modeling of multiple 12 CO, 13 CO 
and C 18 lines by |Dahmen et al. | ( [l998f found that the 
standard conversion factor from the disk overestimated 
molecular gas mass in the Galactic Center by an order- 
of-magnitude. Thus, several independent measurements 
suggest that the Milky Way has aco near its center be- 
tween 3—10 times lower than in the solar neighborhood. 
The properties of the transition from the standard disk 
value to this lower number are not well constrained al- 



though 7-r ays provide some hint of a gradient ( |Strong 
|eTahp004l ). 

The agreement among the various techniques in the 
Galactic Center provides confidence in the resulting low 
aco- On their own, each technique is subject to a num- 
ber of important systematic uncertainties. For 7-ray 
modeling, a key limitation is our knowledge of the cosmic 
ray distribution. In the case of multi-line modeling ap- 
proaches, it is not clear that the approximations involved 
adequately represent the variety of excitation condition s 
in the molecular gas (iMao et al |2000| |Bayet et al.|2006[ ). 
Virial mass based techniques rely on the assumption that 
the clouds are in virial equilibrium, with gravity opposed 
primarily by turbulent motions. If magnetic stresses and 
thermal pressures are significant, or if GMCs are actu- 
ally not self-gravitating, the virial technique may fail to 
reflect the true conversion fa ctor, particularly in r egions 
like the Galactic Center (e.g. Dahmen et al.|[l998 ). 



Measurements of aco with various techniques in 
nearby galaxies have provided somewhat contradictory 
results. Virial mass studies of GMCs in nearby galax- 
ies tend to find aco very similar t o the Milky Way disk 
flWilson|[l995l |Bolatto et al. 2008), even in galaxy cen- 
ters (Donovan Meyer et al. 20 12| and in lo w metallicity 
galaxies like the SMC (Bolatto et al. 2008). In contrast, 
a number of studies using dust- based techniques have 
found very high aco m low- metallicity galaxies, suggest 
ing substantial amounts of "dark molecular gas" or "CO 



free H 2 " around GMCs (Israel 1997; Dobashi et al. 2008 



Leroy et al.||2009b[ |20lip . |Leroy et al. ( |201lJ ) used dust 



to trace the total amount ol gas and determine aco f° r 
the major galaxies of the Local Group. Their results 
suggested that 1) at low metallicity, aco can be well 
above the standard MW value, 2) in the central few kpc 
of M 31, the conversion factor was slightly lower than 
the standard value, and 3) the 10 kpc ring of M 31, 
M 33 and the LMC (which collectively span ~ 0.5 dex in 
metallicity) all showed aco m agreement with the MW 
value. Studies of the centers of several nearby galaxies 
with "large velocity gradient" (LVG) or related modeling 
have found evidence for aco up to an o rder-of-magnitude 



lower than the Milky Way value (e.g. |Weifi et al.j 2001 
Israel||2009b|a 



The study presented here is the first to create maps 
of aco m galaxies outside the Local Group. Therefore, 
we are able to study aco across a range of environments 
and connect the various trends that have been found in 
previous work. We find that for Ico> 1 K km s _1 , aco 
is relatively constant for a range of radii and metallic- 
ities in galaxies. In the central ~kpc, we observe that 
galaxies often show a lower than average aco, sometimes 
by up to factors of 10. These low aco values we find 
in several galaxy centers are not unprecedented, given 
the results from multi-line modeling studies and the ev- 
idence that the Galactic Center shows a similar depres- 
sion. Our lowest meta llicities are still ab ove the transi- 
tion abundance seen by |Leroy et al.| ( |2011[ ) where galaxies 
in the Local Group began to display much higher aco 
(i.e. 12+Log(0/H) < 8.2). The picture that emerges 
from the synthesis of these literature results is one where 
aco is generally around the standard MW value with a 
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factor of ~ 2 variability in the disks of galaxies, in the 
regime where CO is bright and metallicity is greater than 
~ 1/2 Z , with some galaxies showing lower aco m their 
central kpc. 

5.2. Direct Comparison of aco with Virial Mass 
Measurements 

In several cases, our aco measurements overlap with 
previous 12 CO J= (l— 0) virial mass based studies. 



latto et al. 
for 



Bo- 



2008]) derived GMC properties and virial 
UU 2976 and 3077. Their results for 



masses 

N GC 3077 are in goo d a greement with previ ous studies 



by|Meier et al.|fl2001|) a nd Walter e t al.| ( |2002| ). Recently, 
Donovan Meyer et al. ( 2Q12p and Donovan Meyer et al. 
(in prep) have made aco determinations in the central 
regions of NGC 6946 and 4736, respectively, using virial 
masses. In Table 10 we provide a summary of the galax- 
ies in our sample with existing literature measurements. 
Figure 14 shows a comparison between these literature 
virial mass based aco measurements with what is de- 
rived in this study. In almost all cases we find a lower 
aco using the dust-based technique than the virial mass 
technique. 

For NGC 2976, Bolatto et al] ( |2008j ) found conversion 
factors ~ 4 times larger than the Milky Way value on 
average. Almost all of our measurements for NGC 2976 
have aco consistent with the Milky Way value within 
their uncertainties. The central solution pixel has aco 
a factor of 4—5 below the Milky Way value. In con- 
trast, the virial mass and dust based techniques are in 
agreement for NGC 3077, both returning a conversion 
factor consistent with the Milky Way. The center of 
NGC 6946 stands out in our study as places with partic- 
ularly low aco- However, virial mass measurements find 
aco more consistent with the Milky Way value. Fig- 
ure [14] illustrates that the discrepancy decreases with ra- 
dius for NGC 6946 — many of the points with galactocen- 
tric radii closer to r25 ~ 0.2 agree with our dust-based 
measurements. We will discuss the discrepancy between 
the virial and dust-ba sed a co values in the galaxy cen- 
ters further in Section [O) It is interesting to note that 
the normalcy of NGC 3077, currently undergoing a star- 
burst, compared to the also highly star- forming centers 
of 4736 and 6946, suggests that conditions other than 
just high star- formation surface density must contribute 
to determining aco- 

5.3. Direct Comparison with Multi-Line Modeling 
Techniques 

The center of NGC 6946 has been observed extensively 
in molecular gas line s and modeled in a variety of ways. 



Israel fc Baas[(l200l| ) modeled 12 CO and 13 CO lines plus 
[C 1J (492 GHz) observations and found a conversion fac- 
tor ~ 10 times lower than the Milky Way value in the 
center. Using LV G modeling of multiple CO isotopomers 
including 12 C 18 Q, [WaIsh et al.|(|2002[ ) found a conversion 
factor 4 — 5 times below the Milky Way value. Similarly, 
[Meier fc Turner| ( |2004[ ) used an LVG analysis of C 18 



observations, in addition to multiple CO lines, to argue 
that the central region had a conversion factor 4 times 
below the Milky Way value. All of these studies suggest 
a central aco between 4 — 10 times lower than the Milky 
Way value, in good agreement with what we observe with 
dust, but in contrast with the virial mass results. 
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Figure 14. A comparison between virial mass-based measure- 
ments of aco from the literature (shown with red circles) and the 
results of this study. The gray-scale color table for the points from 
this study shows the uncertainty on each aco measurement, as 
in Figure Sample error bars for comparison with the gray-scale 
color table are shown in the top left corner of the plot. Virial mass 
aco measu rements for NGC 2976 and 3077 are taken from |Bo-| 
jlatto et al.|(|2008| >. Measurements fo r NGC 4736 and 6946 are from 
Donovan Meyer et al (in prep) and |Donovan Meyer et al.| ( |2012| , 
respectively. For all of the targets, except JN(JU 3UYY, there is a 
discrepancy between the virial mass and dust based ckco- At larger 
radii, the two measurements are in better agreement in NGC 6946. 

Several other galaxies in our sample have been mod- 
eled with multi-line techniques, however they are at high 
inclination and therefore d o not have reli able aco mea- 
surements from our work. Israel ( 2QQ9b[ ) used multiple 
12 CO and 13 CO lines plus [C I] (492 GHz) observations 
towards the central region of NGC 4631 to constrain aco 
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Table 10 

Literature aco Measurements 



Galaxy 


Regan et al. (2001) 


3cr Below MW 


Literature aco Measurements 




Central Excess? 


oicc)(this study) 


Type 


Reference 


IN^jV^ UDZo 


1ST 
IN 


1ST 
IN 






NGC 2976 




N 


VM 


B08 


NGC 3077 




N 


VM 


B08 


NGC 3351 


Y 


Y 






NGC 3521 


N 


N 






NGC 3627 


Y 


Y 






NGC 4321 


Y 


Y 






NGC 4631 




N 


ML 


109 


NGC 4736 


Y 


Y 


VM 


DM12b 


NGC 5055 


Y 


N a 






NGC 6946 


Y 


Y 


VM, ML 


DM12a, IBOl, MT04, W02 


NGC 7331 


N 


N 


ML 


IB99 


References. 


— B08 - |Bolatto et 


al. (2008); DM12a 


- iDonovan Meyer et al.| (|2012|); DM12b - 



Donovan M eyer et al. (in prep); 109 - Israel (2009b ); IB99"~^|T srael & Baas (1999); IBOl - Israel & 
|Baas| ( 120011); M01 - |Meier et ab] j2UUl| ; MT04 - |Meier & Turner| ( |2U04| ); R01 - |Regan et al.| ^(JETf ; 
W02 - ]Walsh et al.|PPP^ 

Note. — VM - virial mass; ML - multi-line modeling. 

a 2.6cr below MW a C o- 



and found a value 6 ti mes lower than t he standard Milky 
Way value. Similarly, Israel & Baas (|1999| ) argued for 
aco ~ 5 times lower than the Milky Way value in the 
center of NGC 7331. In both of these highly inclined 
galaxies we find aco a factor of 2 or more higher than 
the MW aco- This is due to the failure of our technique 
at high inclinations. 

5.4. Comparison with CO Exponential Disk Profiles 



Regan et al. (2001) found that approximately half of 
their sample of galaxies observed in the BIMA SONG 
survey show excess CO emission in their centers com- 
pared to the extrapolation of their best fit exponential 
disk profile. They argued that this could be due to en- 
hanced reservoirs of molecular gas in the centers or a 
change in the conversion factor. Table [10] lists the cen- 
tral excess classification for th e galaxies in our sample 
studied by |Regan et al. ( 2001[ ). Indeed, we find that in 
every case where they found excess CO above the expo- 
nential profile, we find a conversion factor significantly 
below the Milky Way value. Conversely, for the three 
galaxies we have observed and that they found to not 
have an excess, we find conversion factors consistent with 
the Milky Way value or slightly higher. Similar results 
are found when comparing with the radial profiles de- 
termined from the HERACLES CO maps we utiliz e in 
this study (see radial profiles in Schruba et al.|20lT| ). In 
general, the offset between the central aco and galaxy 
average aco from our measurements is similar in magni- 
tude to the excess observed at the centers compared to 
the extrapolated exponential disk profile. 

6. DISCUSSION 
6.1. Drivers of aco Variations 
In the following, we explore the correlations or lack 
thereof between aco and parameters that may influence 
it. Since our resolution elements are large compared to 
an individual GMC, our aco value is an average over a 
population of GMCs, including any molecular gas that 
may be in a more diffuse phase and not in self-gravitating 
clouds. 



For an individual, self-gravitating, turbulence- 
supported molecular cloud, there are several key 
parameters that influence aco' the density and tem- 
perature of the gas and the fraction of the mass that 
exists in the "CO-dark" phase (where H2 self-shields 
and CO is photodissociated). The general dependence 
of aco o n density and temperature can be illustrated 
using a simple model — a spherical, homogeneous cloud 
where self-gravity is opposed primarily by turbulence 
and 12 CO J =(l— 0) is op tically thick. Following the 
derivation in |Draine (2011), 



aco = 3.4n^(e 



0.5/ p 5.5/T ex 



-1) M pc" 2 (K km s" 



(5) 



Here, 77,3 is the H nucleon density n# = 10 3 n3 cm -3 
and T ex is the excitation temperature. Real molecular 
clouds are certainly not spherical and homogenous, but 
this model allows us to understand the basic dependence 
of aco on density and temperature (note that simulated 
molecular clouds with non-spherica l, non-homogenou s 



structure show similar basic scalings;|Shetty et al. 2011) 



Molecular clouds with increased density will have higher 
aco while clouds with increased temperature have lower 
«co- 

Because CO is optically thick at solar metallicity, the 
abundance of C or O does not directly influence aco- 
However, metallicity may indirectly influence aco by al- 
tering the density or temperature structure of the gas, 
since the heating and cooling rates may be affected by 
metallicity or correlate with metallicity. At lower metal- 
licities, the transition between ionized and neutral car- 
bon to CO may shift relative to the HI/H2 transition 
due to a la ck of dust shielding, leading to layers of "CO- 



dark" H 2 (Tielens fc Hollenbach||1985 
Black|1988 



van Dishoeck 



Wolri re et al.|1993||Kaulman et al.|1999J |Bell 



et al. ]2006f . aco will increase when there are sig nificant 
amounts of gas in these layers. Recent work by |Wolfire 
et al. ( |2010|) used PDR modeling and a spherical model 



lor a turbulent molecular cloud to study variations in 
the conversion factor. They find that the "dark molecu- 
lar gas" fraction (i.e. the fraction of the molecular cloud 
mass where H 2 exists but CO is photodissociated) de- 
pends primarily on the extinction through the cloud or 
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DGR. Studies of simulated molecular clouds have also 
found a dependence of aco on extinction through the 



iouna a dependence 01 acq on € 
cloud ( |Glover fc Mac LowpOlT ). 

The simple model we have described above assumes 
that the cloud's self-gravity is balanced by turbulence. 
If there are other forces that contribute to the virial 
equation, aco m ay be altered. For a self-gravitating 
cloud, an increase in the pressure on the surface of the 
cloud will lead to a larger velocity dispersion. Because 
CO is optically thick, the amount of emission that "es- 
capes" in the CO line is directly tied to the velocity 
dispersion of the cloud and increasing the velocity dis- 
persion increases the CO emission for the same mass of 
H 2 . Therefore, if a self-gravitating cloud were subject 
to significant external pressure we would expect a lower 
q;co- Conversely, if magnetic fields played an important 
role in supporting the cloud, this would lower the veloc- 
ity dispersion and raise ftco- In ultra-luminous infrared 
galaxies (ULIRGs), the molecular gas may not be in indi- 
vidual self-gravitating clouds, but instead in a large-scale 
molecular medium where stellar mass contributes to the 



diation field may play a role in determining T ex ( |Womre 
et al.|[l993 ). Likewise, qpAH and metallicity could both 



gravitational potential, (e.g. |Downes fc Solomon 
In this case, the velocity dispersion is larger than what 
would be expected from the gas mass alone and the CO 
emission is enhanced for a given gas mass (resulting in a 
lower aco)- 

In a ~kpc region of a galaxy, many individual molecu- 
lar clouds will be averaged together in our measurement. 
Therefore, another factor which could contribute to the 
variation of our measured aco is changes in the cloud 
populations. In addition, any significant component of 
diffuse CO emission would be included as well. Several 
studies have suggested that CO emission from diffuse gas 
may not be negligible, even in the local area of th e Milky 
Way ( [Liszt fc Lucas|1998| [Goldsmith et a"Ll|2QQ8| ). Inter- 
estingly, it appears that lor the Milky Way, aco 1S si milar 



in diffuse molecular gas and self-gravitating GMCs (Liszt 
et al.||2010 ). However, it is possible that under differing 



environmental conditions the conversion factors appro- 
priate for these two types of gas could be quite different. 
The contribution of diffuse molecular gas with a different 
aco has been cited previously as a cause for discrepan- 
cies between different techniques for measuring aco in 
M 51 ( |Schinnerer et al.||2QlQ[ ). 

To summarize our theoretical expectations, aco mea- 
sured in ~kpc regions of nearby galaxies can vary as 
a function of environment for many reasons: changes 
in the excitation temperature and density of the gas or 
contributions from pressure or magnetic support in self- 
gravitating clouds; envelopes of "CO-dark" gas; changes 
in the molecular cloud population; or contributions from 
diffuse CO emission. Unfortunately, directly measur- 
ing temperature, density and velocity dispersion in large 
samples of extragalactic GMCs is challenging due to the 
need to resolve individual clouds in multiple molecular 
gas emission lines. Therefore, we are left with more indi- 
rect tracer s of potential changes to the GMC properties. 
In Section [4~6] we examined the correlation of aco with 
U, c[pah, metallicity, £*, Hsfr, £d and galactocentric 
radius. 

The average radiation field U could influence aco 
through the gas excitation temperature. If photoelectric 
heating dominates over other heat sources (e.g. cosmic 
rays) at the tco ~ 1 surface, then the intensity of the ra- 



potent iaTIy influence the efficiency of the pho toelectric 
effect flBakes fc Tielens|1994l iRollig et al.|2 QQ6), thereby 
changing the heating rate. We expect that £sfr should 
be responsible for higher U in many regions. Enhanced 
SFR leading to higher radiation field intensities has been 
suggested as the explanation for high obse rved aco in 
sever al outer-disk molecular clouds in M 33 ( |Bigiel et al. 
2010|). If the gas excitation temperature were affected 



by the radiation field or photoelectric heating, we would 
expect negative correlations (i.e. lower aco at higher U, 
^sfr, qpAH, etc). The correlations we observe between 
aco and ?7, qpAH and Ssfr are generally weak. U and 
aco show the strongest association and the slope of the 
trend is negative, with lower aco at higher U. This is 
potentially consistent with the expectation of higher ra- 
diation field intensities leading to warmer molecular gas 
temperatures, but the correlation is weak and other vari- 
ables may play a more important role. 

We see a weak trend of aco with metallicity. Since 
our observations are limited to regions with metallicities 
similar to or higher than the Milky Way, we may not ex- 
pect to see a strong correlation between aco and metal- 
licity. The fraction of gas mass in the layer where CO 
is photodissociated is p redicted to be ~ 3 0% at Milky 
Way metallicity/DGR flWolfire et aL]|201Q| ). Thus, in- 
creasing the DGR should only have a minimal effect on 
the conversion factor at MW metallicity and above. To 
illustrate this point, in Figure [15| we plot our aco mea- 
surements a s a fun ction of metallicity and overlay the 
|Leroy et a l.l (|201l | Loca l Grou p me asurements and the 
model s ol | Woltire et al. ( 2Q1Q[ ) and Glover & Mac Low 
( 201 1 ) . The model predictions are normalized to have 
a C o= 4.4 M pc" 2 (K km s" 1 )" 1 at the MW metallic- 
ity (i.e. 12+Log(0/H) = 8.5 in PT05 and 8.8 in KK04; 
Croxall et al., in prep). Both models assume a linear scal- 
ing of DGR with metallicity and a fixed gas mass surface 
density for molecular clouds. The model predictions are 
discussed in detail in Leroy et al. (submitted). In both 
calibrations, our measurements do not extend into the 
regime where "CO-dark" H2 dominates aco- We note 
that an investigation of the aco-DGR relationship using 
our results would require taking into account the strong 
correlation of the uncertainties in these parameters. 

The strongest correlation we observe for aco is that 
with £*. This correlation is the only one that is stronger 
than the correlation of aco with galactocentric radius 
T25. In essence, our definition of T25 is based on a scale 
length of the stellar disk, defined by a B-band surface 
brightness, so it is possible the stronger correlation of 
aco with is due to it being a better proxy for the 
disk scale length than 125. In addition, of the environ- 
mental parameters we have available for our analysis, £* 
is the most straightforward to measure and may show 
the smallest systematic uncertainties. It may be the 
strongest correlation simply because the other parame- 
ters are more uncertain. The question remains, however, 
as to why aco would correlate with the stellar mass sur- 
face density. There are a number of possibilities: pres- 
sure contributions to GMC virial balance, enhanced frac- 
tions of diffuse CO emission correlated with ISM pres- 
sure, changes to the population of GMCs, or other effects. 
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Figure 15. aco as a function of metallicity compared to previous measurements and models. The left panel shows measurements in 
the PT05 calibration while the right panel shows KK04. Measurements from this paper are shown with gray points (individual solutions) 
and red circles (galaxy averages). We show all of the solution pixels where Ico> 1 km s _1 and i > 65°, regardless of the source of the 
metallicity measurement (i. e. without the requ irement of having measured HII region metallicities from M10). Measurements of aco m 
Local Group galaxies from | Leroy et al.| ( |2011| are shown with green circles. The MW aco is shown with a g ray line and the av erage 
of our solution pixels with no weighting is shown with a dotted black line. Predic tions based on the model of Wol fire et al.| (2010) are 
shown with a purple dot-dashed line and those based on |Glover Mac Low| ( |20ll| ) are shown with a dashed blue line. These predictions 
assume a linear dependence of DGR on metallicity and a fixed gas mass surface density for molecular clouds of Sqmc — 100 Mq pc -2 . 
The model predictions are normalized to have aico= 4.4 Mq pc -2 (K km s -1 ) -1 at the metallicity adopted for the Milky Way in each 
calibration. Note that the metallicities we use for NGC 5457 are not from strong line calibrations, so it appears in the same position both 
plots. Regardless of the metallicity calibration, our measurements do not extend to low enough metallicities to constrain the effects of 
"CO-dark" H 2 . 

We do not speculate here about what causes this corre- 
lation, only note that observations of multiple molecular 
gas lines at GMC resolution will help understand these 
trends and are possible with ALMA. 

Although we see some evidence for weak correlations 
between aco and environmental variables, our primary 
result is that in regions with Ico> 1 km s _1 in these 
galaxies, aco is mostly insensitive to environment except 
in the central regions. This may be due to our limited 
S/N for each individual aco measurement plus the rela- 
tively small range of environmental conditions we probe. 
Another possibility is that in the disks of most galaxies, 
the properties of molecular gas are insensitive to envi- 
ronment because most of the gas resides in molecular 
clouds that are not under significant external pressure or 
subject to enhanced turbulence or magnetic fields. The 
conversion factor appropriate for such clouds should be 
similar to the standard Milky Way value. 



6.2. Comparison with Galaxy Simulations 

Because of the resolution of our observations, our mea- 
surements average over a population of molecular clouds. 
This means that environmental variations in the cloud 
population can also be responsible for changes in aco- 
Therefore, it is interesting to compare our results to what 
has been found with galaxy simulations, despite the fact 
that these simulations are forced to adopt a sub-grid 
model that prescribes aco for unresolved clouds. In Fig- 
ure [16] we show our measured aco values as a function of 
the average Ico in each solution pixel. Due to our com- 
pleteness cut, we have no points below 1 K km s _1 . The 
mean and standard deviation of the aco measurements 
in bins of 0.25 K km s _1 are shown with black circles and 
error bars. 

Overlaid on Figure [16] we show the p redict ions of two 
recent simulations. |Feldmann et al. ( |2Q12| F12) cou- 



ing a fixed gas temperature and either a c onstant line 
width of 3 k m s _1 or a virial scaling. The Narayanan 
et al. ( 2Q12| N12) predictions include sub-grid semi- 
analytic models for the chemical and thermal state of 
the gas and dust. In the N12 models, the temperature 
of the molecular gas is not fixed and depends on heating 
by cosmic rays and, in dense regions, on energy transfer 
with dust. For both models it is necessary to adopt a 
filling factor or clumping correction due to the difference 
in resolution between the simulations and our observa- 
tions. F12 and N12 both have ~ 60 — 70pc resolution, 
so we apply the same correction to both. For purposes 
of comparing with the trend of our observed aco as a 
function of CO integrated intensity, we have adopted a 
filling factor correction of 3. This correction factor is 
larger than purely the ratio of the resolutions because 
our beam will contain more than one molecular cloud. 
Verifying that this absolute scaling is correct would re- 
quire a more detailed investigation that is beyond the 
scope of this paper. 

We find that the predicted trend of aco wit h CO in- 
tegrated intensity from |Narayanan et al.| (|2012) matches 
the average observed behavior well, aside from under- 
predicting aco in the central regions of some galaxies. 
There is c onsiderable scatt e r aro und this trend, how- 
ever. The jFeldmann et al.| p012) results diverge from 
our measured trend at low (JO integrated inten sities, pre- 
dicting a ste eper increase than we observe. Feldmann 
2Q12[ ) argue that on ~kpc scales metallicity is 
primary driver of the conversion factor variations. 
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Figure 16. A comparison bet ween our ar,r> me asurement s and 
the predictions of si mulations bylNarayanan et al.|(2012) and Feld-| 
|mann et al.| (|2012| >. On the x-axis we show the average Icq m 
a given solution pixel. All solutions above our Ico an d inclina- 
tion cuts are shown as gray circles. The simulation predictions 
are shown at three possible metallicities: Z = Zq in a solid line 
and a factor of 2 above and below that with dashed lines. The 
mean and standard deviation in 0.25 K km s — 1 bins are shown 
with black circles and error bars. The predictions from the simu- 
lations have been divided by a filling factor of 3 to correct for our 
different resolutions. In general, the average behavior of our ago 
measuremen ts agrees well with the predictions of the Narayanan 
et al.|(|2012| m odel except at the highest Ico- 



mann et al. (2012) assume a metallicity dependent DGR. 
Our aco results do not show a strong dependence on 
metallicity, in contrast to this prediction. One key issue 
may be the assumption of a fixe d gas temperature in the 
clouds. Narayanan et al. (2012) do not make such an as- 
sumption and see a decrease m aco in regions with high 
star-formation rate surface densities due to enhanced gas 
temperatures and velocity dispersions. 

6.3. Low Central aco an d Discrepancies with Virial 
Mass Based Measurements 

For some galaxy centers in our sample we measure aco 
lower than the typical Milky Way value by factors of 
5—10. Interestingly, the virial mass based aco measure- 
ments for several of these regions do not agree with the 
values found here, instead recovering aco values close 
to that of the Milky Way. Results of multi-line mod- 
eling are generally in much better agreement with our 
results. These centers also tend to be regions where the 
C O shows an excess o ver the exponential disk, as found 
by Regan et al. (2001 ). From Figure [l4| we found that the 
discrepancy between virial mass and dust-based aco de- 
creases with radius for NGC 6946 — suggesting that this 
is a property of the galaxy centers rather than an issue 
with the resolution of the virial mass techniques (which 
tends to be comparable to the size of the GMCs). 

There are several scenarios in which aco m &y decrease 
but the virial mass based aco would appear similar to 
the MW value. If the molecular gas is in self- gravitating 
GMCs, it is possible to decrease aco by increasing the 
velocity dispersion of the cloud due to the effect of exter- 
nal pressure. In such a situation, the measured line- width 
would lead to a larger virial mass and the cloud's CO 
luminosity would also increase. Thus, by comparing the 
virial mass and CO luminosity of more highly pressurized 
clouds, one could derive a conversion factor similar to the 



MW value, driven by the fact t hat the virial mas s has 
overestimated the true H 2 mass. Oka et al. (1998a) have 
suggested that external pressure in the Galactic Center 
can account for up to an order-o f-magnitude change in 



the c onversion factor. However, Donovan Meyer et al. 
( [2012] ) found that GMCs in NGC 6946 did not gener- 
ally show enhanced velocity dispersi on, in contrast to th e 
Galactic Center clouds identified by Oka et al. ( 1998b[ ). 

Another potential explanation for this discrepancy is 
that gas in the centers of some galaxies is not in self- 
gravitating GMCs, but rather in a larger-scale molecular 
medium that is gravitationally bound to the stars and 
gas. In this situation, the velocity dispersion is larger 
than would be the case for self-gravity of the gas alone, 
and more CO emission is produced for a given amount 
of gas . This expl anation has been suggested by |Downes| 
|et al.| ( |1993| and Downes fc Solomon| ( |1998[ ) to explain 
the properties of molecular gas in the nuclear disks of 
ultra-lumino us infrared galaxies (ULIR Gs) (see also re- 
cent work by |Papadopoulos et al.| 2012). 

These scenarios should alter the velocity dispersion, 
temperature and optical depth of the molecular gas. 
With the resolution of the available data, we cannot ob- 
tain reliable measurements of velocity dispersion in the 
galaxy centers. Each resolution element encompasses a 
range of velocities purely due to rotation, so we cannot 
disentangle dispersion from rotation. The resolution is 
also far too low to isolate individual GMCs. We do have 
some hope of testing predictions about the optical depth 
and temperature by examining the ratios of 12 CO and 
13 CO lines. Usero et al. (in prep) presents measurements 
of 13 CO and 12 CO lines towards regions of several HER- 
ACLES galaxies. Selecting all of their pointings towards 
the disks of galaxies in our sample (in regions above our 
Ico and inclination cuts), the uncertainty weighted mean 
of #12/13 = 12 CO J=(1-0)/ 13 CO J=(l-0) is 8.24±0.11. 
The only galaxy centers with 13 CO measurements are 
NGC 0628, 3184, 5055 and 6946. NGC 0628, 3184 and 
5055 show #12/13 consistent with or slightly lower than 
the average disk value. NGC 6946, on the other hand, 
has #12/13 = 18.53 ±0.81. The high #12/13 for the center 
of NGC 6946 has been noted in several previous studies 
flPaglione et al.|2001| [Israel fc Baas|2001| Meier fc Turner| 
|2004[ ). Ol those lour galaxies, NGC 6946 is the only one 
with a central aco more than 3a below the MW value. 
Further attempts to explain the low aco values in some 
galaxy centers will require more extensive mapping in 
other molecular gas lines, preferably at high spatial res- 
olution. 

6.4. Dust-to-Gas Ratio 



As shown in Figures 12 and[13j we observe a good cor- 
relation of DGR with metallicity. Unlike most previous 
studies of DGR resolved within galaxies, we made no 
assumption about aco m determining these values. Al- 
though we do not probe a wide range of metallicit y d ue 



to the limitations on CO S/N, it is clear in Figure 13 



particular, that DGR is correlated with metallicity with 
less than a factor of 2 scatter over 0.5 dex in metallic- 
ity. A linear dependence of DGR on metallicity suggests 
a constant fraction of heavy elements are locked up in 
dust grains. Chemical evolution and dust life-cycle mod- 
els have varying predictions for the dependence of DGR 
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on metallicity jDw ek 1998 Lisenfeld fc Ferraraj 
rashita et al. 2002p . Most of the dramatic cfi 



Ferrara|1998|p^ 
dramatic differences 
between the models, however, occur at lower metallici- 
ties and observations of dwarf galaxies may pro vide more 



leverage on d istinguishing between models (e.g. 
|Camus et al1|2Q12| . 



Herrera- 



6.5. Recommendations on Choosing aco 

Due to the requirements on CO S/N to apply the tech- 
nique we have used, our measurements only tell us about 
aco above Ico=l K km s _1 in the galaxies we have 
targeted. Generally, this limits us to the inner parts of 
galaxies (<r2s) where metallicities are comparable to the 
MW disk. However, these are the regions where H2 con- 



tributes m ost significantly to the total gas mass ( Schruba 
|et al.pQlTj ), so applying our conversion factor results to 
unresolved galaxies should work reasonably well for re- 
covering the total H2 mass. A forthcoming paper ex- 
tending our study to the HI dominated regions of the 
KINGFISH galaxies, will test this assertion. Note that 
our recommendations are only applicable to spiral galax- 
ies with metallicity similar to that of the MW. 

For an unresolved galaxy, we recommend adopting our 
galaxy-based average aco of 3.1 M pc -2 (K km s -1 ) -1 
with uncertainty of 0.3 dex (a factor of ~ 2). This is the 
mean and standard deviation of the average value for all 
galaxies with inclination less than 65°. 

In dealing with resolved galaxies, we recommend 
adopting a flat radial aco profile for regions with Ico 
> 1 K km s _1 except in the central ~ 0.1r25, where the 
average aco is a factor of 2 lower. When dealing with 
a single galaxy, our average radial profile from Figure [7] 
suggest that the aco values have ~ 0.2 dex of scatter in 
each radial bin (0.3 dex for the central bin). Therefore, 
for studies of molecular gas within a single galaxy, the 
relative values of aco adopted from this profile have a 
factor of ~ 1.5 uncertainty. When comparing resolved 
galaxies, it is necessary to further recognize that the ab- 
solute normalization of the radial profiles, i.e. the galaxy 
average value, has an additional 0.3 dex of uncertainty. 

In the case where additional information about the 
galaxy is available, such as profiles of stellar mass sur- 
face density, star formation rate surface density or line 
ratios of 12 CO (2 - 1)/(1 - 0) or 12 CO/ 13 CO, we sug- 
gest using these maps to pick out regions that may have 
aco ver y different from the mean. This includes regions 
with > 1000 M pc" 2 , E SF R > 0.1 M yr" 1 kpc" 2 , 
or R21 ^ 1 at our working resolution of ~ 1 kpc. In 
these regions, aco can be 5—10 times lower than the 
MW aco, and are often systematically lower than the 
average radial profile. 

7. SUMMARY & CONCLUSIONS 

The availability of high angular resolution far-infrared 
maps from the Herschel Space Observatory and sensi- 
tive 12 CO J=(2— 1) and HI maps has recently allowed us 
to trace dust and gas mass surface densities in nearby 
galaxies on ~kpc scales. On these scales, we expect (and 
have verified) that the DGR is approximately constant. 
Therefore, we are able to combine the dust mass sur- 
face density maps with matched resolution CO and HI 
to solve simultaneously for aco and DGR. The solution 
technique finds the aco that best minimizes the scat- 
ter in DGR values in ~kpc regions across the galaxies. 



We have performed a thorough investigation into the ef- 
ficacy of this technique, using data from the KINGFISH 
key program on Herschel, the large IRAM 30m survey 
HERACLES and the THINGS HI survey with the VLA. 
Our tests show that above Ico ~ 1 K km s _1 we reliably 
achieve accurate solutions for aco- We have used a fixed 
ratio between the (2—1) and (1—0) lines to present our 
aco results on the more typically used 12 CO J=(l— 0) 
scale and have shown using literature measurements that 
variations in the line ratio do not effect our results in 
these galaxies. 

We find that the average aco fc> r the galaxies in our 
sample is 3.1 M pc -2 (K km s -1 ) -1 . This value is 
slightly lower than the Milky Way aco = 4.4 M pc -2 
(K km s -1 ) -1 , but the MW value is well within the stan- 
dard deviation of our measurements. Treating all 782 
solutions independently, we find aco = 2.6 M pc -2 (K 
km s -1 ) -1 with 0.4 dex standard deviation. In all aver- 
ages, we have removed galaxies with inclinations higher 
than 65° since our solution pixel samples gas along the 
line-of-sight with a range of DGR and therefore does not 
conform to our main assumption — that the DGR is con- 
stant in the solution pixels. 

Within the galaxies we observe a relatively flat profile 
of aco as a function of galactocentric radius aside from 
in the galaxy centers. Normalizing each galaxy by its 
average aco the average galaxy shows a factor of ~ 2 
lower aco in its center. In several galaxies, this central 
value can be factors of 3 or more lower than the galaxy 
average. In several notable cases, the central aco value is 
factors of 5 to 10 lower than the standard MW aco= 4.4 
M pc" 2 (K km s- 1 )- 1 . 

We have investigated the correlations between aco and 
environmental parameters in an attempt to isolate fac- 
tors that drive variations in aco- The strongest corre- 
lation we find is between aco and stellar mass surface 
density (£*). If the strength of this correlation relative 
to the other parameters is real and not due to issues 
with measuring the other parameters, there are a num- 
ber of potential explanations for why aco should depend 
on stellar mass surface density, including the influence 
of ISM pressure on both molecular clouds and a more 
diffuse molecular medium. Distinguishing among these 
possibilities will require ALMA observations that resolve 
individual GMCs in multiple molecular gas lines across 
a range of extragalactic environments. 

We do not observe a strong correlation between aco 
and metallicity in our galaxies. In the range of metallici- 
ties we have sampled, this conclusion is not unexpected. 
The abundance of C or O does not play a major role in 
determining aco at these metallicities since CO is opti- 
cally thick in molecular clouds. The metallicity primar- 
ily influences aco through the effects of dust shielding 
(due to the dependence of DGR on metallicity). Mod- 
els and simulations have suggested that dust shielding 
can influence the amount of gas in "CO-dark" layers of 
GMCs, but at MW metallicity, only ~ 30% of the gas is 
in this layer. At lower metallicities, a strong dependence 
of aco on metallicity has been observed using the same 
technique we employ ( Leroy et al.||2011 ). 

In the centers of several galaxies, we find aco values 
5—10 times lower than the MW aco- These regions are 
also significantly below the average aco f° r their galaxy. 
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Comparison of our measured aco with values from the 
literature shows good agreement between our dust-based 
results and multi-line modeling results. In contrast, our 
aco values in these regions can be much lower than aco 
determined from virial mass based techniques. The dis- 
crepancy with virial mass measurements becomes smaller 
at larger galactocentric radii, suggesting this is a partic- 
ular property of the gas in galaxy centers. At the res- 
olution of our observations, the central region with low 
aco is unresolved (< kpc). We suggest several potential 
explanations for the low aco value and the fact that it is 
not reflected in virial masses including ISM pressure con- 
tributions to GMC virial balance, increases in molecular 
gas temperature and/or a more diffuse molecular medium 
similar to what is found in ULIRGs. With the limited 
amount of 13 CO and other molecular line observations 
in this region, it is not possible to distinguish between 
these scenarios, but NGC 6946 at least shows some evi- 
dence for potentially lower CO optical depth towards its 
center. As is explored in Leroy et al. (submitted), these 
galaxy centers show enhanced star-formation efficiencies 
when we apply our aco, close to what is seen along the 
"sta rburst sequence" in ultra-luminous infrared galaxies 
(e.g. |Daddi et al.|201Q[ ). 

In addition to aco we also simultaneously measure 
DGR in all of our regions. On average, we find 
Log(DGR)= —1.86 treating all solution pixels equally, 
with a standard deviation of 0.22 dex. When we force 
each galaxy to contribute equally, the average is es- 
sentially indistinguishable from the MW DGR. Unlike 
aco, DGR is well-correlated with metallicity, with a 
slope slightly shallower than linear (although this slope 
is mainly due to NGC 3184 being offset from the main 
trend, removing it from the sample gives a slope consis- 
tent with linear). The approximately linear dependence 
of DGR on metallicity agrees with the predictions of dust 
evolution models, but our measurements do not cover 
a wide enough metallicity range to distinguish between 
them. 

The results presented here suggest a picture where 
aco is slightly lower than the typical value of 4.4 M 
pc -2 (K km s -1 ) -1 in the disks of most galaxies, and 
mainly constant as a function of radius despite changes 
in metallicity, radiation field intensity and star-formation 
rate surface density. Galaxy centers appear to be a 
different regime, where external pressure or changes in 
the character of molecular gas (i.e. mostly confined 
to self-gravitating GMCs versus a more diffuse molec- 
ular medium) may bring about large changes in aco- 
Through the galaxy, however, DGR appears to be an 
approximately linear function of metallicity. The simple 
behavior of DGR provides a unique tool to study the 
ISM in nearby galaxies, if we can obtain measurements 
of metallicity gradients with trustworthy calibration. 
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APPENDIX 

TECHNIQUE FOR SIMULTANEOUSLY DETERMINING DGR AND a C o 

Prior to choosing the technique described in Section [3] we have explored a variety of possible ways to simultaneously 
constrain DGR and aco- The various techniques represent, in essence, different ways to judge how uniform the DGR 
is in the region in question, at any given value of aco- After defining these variations on the essential technique, we 
construct simulated datasets with known aco and DGR and realistic statistical errors and intrinsic scatter to test 
the solution techniques. These simulated data tests allow us to optimize the performance of the technique given the 
properties of our dataset. 

Potential Techniques 

The general procedure for finding the most "uniform" DGR in the solution pixel given the spatially resolved mea- 
surements of Shi, and Ico is to step through a grid of aco values, calculate DGR at each sampling point in 
the solution pixel using that aco, measure some statistic that can be used to judge how uniform the DGR values 
are across the region and finally, either minimize or maximize that statistic as a function of aco- Below we describe 
the various statistics we have tested to judge the uniformity of the DGR values in a given region. The reason we 
investigate all of these possible statistics is because each weights outliers differently, and it is not obvious a priori 
which will work best to recover the underlying aco- For techniques that involve minimizing scatter in the DGR, we 
use three different statistics to measure scatter: the standard deviation (RMS), t he median absolut e deviation (MAD) 
and a robust estimator of the standard deviation from Tukey's biweight (BIW) ( |Press et al.||2QQ2[ ). 

1. Minimum Fractional Scatter in DGR (FS): At each value of aco we calculate the DGR for all points in the 
region in question. We then divide each DGR value by the mean DGR in the region and measure the scatter of 
the resulting values (with RMS, MAD and BIW statistics). 

2. Minimum Logarithmic Scatter in DGR (LS): At each aco we calculate the scatter of the logarithm of the DGR 
values (with RMS, MAD and BIW). This is similar to the minimization of fractional scatter, but with a different 
weighting for the individual points. This also throws away data points with negative values, since their logarithm 
is undefined. 

3. Maximum Correlation Coefficient of Ed and Ec a s (LC or RC): At each value of aco we compute the linear and 
rank correlation coefficients between Ed and E gas = Ehi + acolco- The best aco value will maximize the 
correlation between these two quantities. 

4. Minimum x 2 of Best-Fit Plane to Ico, Nhi and Ed (PF): The Ed, Ehi and Eh 2 values describe a plane with 
two free parameters. At each value of aco we fit a plane to the measurements to determine DGR. The best 
value of aco will De the one that minimizes the x 2 of the best-fit plane. 

5. Minimum Correlation Coefficient of DGR vs. Ico/Nhi (CHC): In this technique we search for the aco value 
that minimizes the correlation of DGR with Ico/Nhi- This technique weights the points differently depending 
on how much local E gas depends on aco- 

Definition of the "Solution Pixel" Size and Location 

A key aspect of the technique we use to constrain DGR and aco is defining the region over which we assume there 
to be one value each of DGR. Under ideal circumstances, the regions would be as small as possible given the resolution 
of the maps. Realistically, however, there are several considerations that must be taken into account when defining the 
regions including: covering an adequate range of CO/HI ratios, having enough points to measure scatter accurately, 
and keeping the region to ~kpc scales over which it is reasonable to assume that the DGR is indeed constant.. 

To complement the hexagonal, half-beam spaced grid with which we sampled the original data (i.e. our "sampling 
points") we again use hexagonal spacings to define the "solution pixels" as the individual samples can be divided 
naturally into concentric hexagons. We tested a variety of solution pixel sizes ranging from 19-point to 271-point 
hexagons. 

We found that we 37-point solution pixels was a good compromise between size and solution quality. Because of 
the underlying half-beam spaced sampling grid, the 37 point solution pixel contains ~ 9 independent measurements. 
For the next smallest possible hexagon (19-points), there are not sufficient numbers of independent measurements to 
reliably judge the scatter in the DGR. 

In order to fully sample the map, the center of each hexagonal pixel is offset by 1/2 of the spacing from its neighbor. 
This results in the "solution pixels" not being independent-each neighboring pixel shares ~ 40% of the same data (the 
oversampling can be exactly described as n 2 /(3n 2 — 3n + 1) for concentric hexagons). 
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Test of the Solution Technique with Simulated Data 

We have investigated the efficacy of the various techniques using simulated data modeled on our observations. Our 
goal is to identify the technique or techniques that return the most accurate values of ftco and DGR over the range 
of conditions and signal-to- noise levels in our dataset. We also aim to characterize any biases in the recovered ftco as 
a function of input aco and S/N. In performing these tests, it is important to model the simulated data on our real 
observations. This is because the ability of the techniques to judge the DGR uniformity depend not only on the S/N 
of the measurements but also on the range of CO/HI ratios. Since CO and HI can be correlated, anti-correlated or 
independent of each other on kiloparsec scales, depending on the state of the ISM, it is difficult to generate a reasonable 
set of simulated data from scratch that encompasses the range of CO/HI behaviors in the observations. Therefore, we 
expect the best test of the technique will come from simulated data that are closely modeled on the observations. In 
the following, we describe how the simulated data are generated and discuss the results of the test. 

The simulated data are generated via a Monte Carlo procedure in which we randomly choose a galaxy from our 
sample (listed in Table [I]) at SPIRE 350 resolution, choose a solution pixel and create simulated data based on the 
selected CO, HI and Ed samples. For each trial we randomly sample a range of aco values between 0.5 — 50 M 
pc -2 (K km s -1 ) -1 and intrinsic scatter per axis (evenly divided among the Ico? ^hi and £d axes, for lack of better 
knowledge of the sources of intrinsic scatter) between 0.0 — 0.4 dex. We note that because of the properties of lognormal 
distributions, the scatter in the DGR is not simply V^x the input intrinsic scatter per axis, so we approximate it 
after-the-fact using the input scatter per axis and the mean H2/HI ratio for the simulated measurements. 

Our procedure is the following: 1) select the observations with the sampling described above, 2) use the Ico 
and Shi measurements along with the random aco value to create a simulated S gas vector, 3) calculate the DGR 
value that preserves the average S/N level of the Ed measurements and create a simulated, noise- free Hd vector by 
multiplying the S gas from step 2 by this DGR, 4) apply the randomly-selected intrinsic scatter to each axis and 5) add 
simulated measurement errors according to the observed uncertainties in the Ico? ^hi and £d maps. These simulated 
measurements are then used as input to solve for aco and DGR using the techniques described above. 

From this Monte Carlo simulation, we can then choose the technique that is the most robust (i.e. most frequently 
achieves a solution) and accurate (i.e. to get closest to the known input aco)- In Figure [TF| we show the results of 
this simulation for the techniques we have outlined. The panels show (left) the mean difference between the input and 
recovered aco value, (middle) the standard deviation of the difference of the input and recovered aco, and (right) the 
percent of the Monte Carlo trials that achieve a solution. These values are binned in a two-dimensional space defined 
by the intrinsic scatter added to the simulated data and the median S/N of the simulated CO measurements. 

It is clear that the FS and LS techniques provide the least bias in the recovered aco over the relevant range of intrinsic 
scatter and S/N (except for the FS_RMS which is significantly biased at low CO S/N). The correlation coefficient (LC 
and RC), planefit (PF) and CO/HI correlation (CHC) techniques all provide much more limited regions with low 
bias in the recovered aco as well as having a lower fraction of solutions. In general the MAD based techniques have 
slightly less bias towards the low CO S/N region, but show a larger scatter between the input and recovered aco- 
Overall, we find that the LS_BIW (robust biweight mean of the logarithmic scatter) is the technique that shows the 
best performance in accuracy, precision and robustness over the range of characteristics of our dataset. 

In addition to helping us choose the best technique, the tests with simulated data also allow us to investigate how 
the measured minimum in the logarithmic scatter corresponds to the known intrinsic scatter added to the simulated 
data. In addition to being an interesting quantity scientifically, the measured scatter can also help us to identify which 
region of the parameter space a given solution falls in, allowing us to judge how biased we expect it to be given our 
knowledge from these simulated data tests. Figure [TB] shows a comparison of the known input intrinsic scatter with the 
measured minimum in the logarithmic scatter for the LS_BIW technique. Most of the outliers come from regions with 
low CO S/N. At very low levels of intrinsic scatter, the measured scatter is larger due to the influence of observational 
uncertainties. In general, however, the measured minimum log scatter is a very good proxy for the intrinsic scatter in 
the DGR. 

Based on these results with simulated data, we can use the measured median CO S/N and the measured minimum of 
the scatter in the DGR values to predict how biased and how uncertain a given determination of aco f° r °ur dataset. 
To do this we use the measured minimum in the DGR scatter as an estimate of the intrinsic DGR scatter (it is biased 
slightly higher, but this is a minor effect and represents a conservative estimate of the intrinsic scatter). Then using 
this value and the median CO S/N in the solution pixel we interpolate in the grid presented in Figure [Tt] for "LS_BIW" 
to predict the bias inherent in the technique. 

Testing for Potential Bias at High CO S/N and Low Shi 

One of th e key results of our work is the measurement of low aco values in the centers of some galaxies (see 
Section [ob. These locations are outliers from the average solution pixel in the sense of having high CO S/N and often 
a relative deficiency in Hhi- Since they make up a small subset of all solution pixels, our Monte Carlo simulation from 
the previous Section may not adequately judge the bias of the technique in these regions, since they are not sampled 
frequently enough (only 26 out of ~ 900 solution pixels). We therefore perform a separate Monte Carlo simulation 
for the galaxy centers in order to verify that these low aco values are not due to biases in the technique (we have 
demonstrated that the average pixel in the sample is not biased using the experiment in the previous section). 

The results of the technique are biased if it returns an aco systematically different from the true aco- For this 
test, we will define several different aco values. First, aco,meas is the measured value for a given solution pixel, we 
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Figure 17. Results of simulated data tests to optimize the solution technique. Columns correspond to: left, the difference between 
the known input aco and the recovered best-solution aco (A«co); middle, the standard deviation of the difference between input and 
recovered aco (o - ); an d right, the fraction of trials in each bin where a solution is achieved. The results are presented as two-dimensional 
histograms, binned in CO S/N and the simulated intrinsic scatter added to the data. The rows show results for the various solution 
techniques described previously. From top to bottom in the left panel these are: minimization of the fractional scatter using the mean 
(FS_RMS), median (FS_MAD) and robust mean (FS_BIW) and the minimization of the logarithmic scatter using the mean (LS_RMS) 
and median (LS_MAD). In the right panel we show minimization of the logarithmic scatter using the robust mean (LS_BIW), minimizing 
the linear correlation coefficient (LC), minimizing the rank correlation coefficient (RC), minimizing the % 2 of the best-fit plane (PF) and 
minimizing the correlation coefficient between DGR and Ico/^HI (CHC). The contours in the Aaco column show the regions of parameter 
space where the difference in the input and recovered aco is l ess than 0.2 dex. The LS_BIW technique minimizes the difference between 
input and recovered aco an d achieves a small a over the largest area, although the distinction from FS_BIW is not large. 

do not know if this measured value is biased, it is only what we have measured for that pixel given the true data. 
Second, aco, in is a known aco value we have used to generate simulated data. Finally, aco, out is the aco value that 
we recover for aco, in using our solution technique. 

The essence of this test is that we take a range of aco, in values and map them to their aco, out values, by generating 
simulated data given aco, in and running it through our solution technique. This simulated data generation proceeds 
as described in the previous section and is designed to preserve the CO/HI ratio and S/N of the measurements for a 
given pixel. In addition to observational noise, we also add intrinsic scatter between 0.01—0.1 dex to the simulated 
data, which reflects the level of intrinsic scatter we have found in the real solutions. For each solution pixel we generate 
10 4 random aco, in values that span our range of aco= 0.1 — 100. 

With this simulation, we can determine for any given value of aco, out? which values of aco, in ended up there. If a 
solution is biased, the correspondence between aco, in and aco, out will not be one-to-one. For example, a bias may 
make us consistently recover aco, out = 1 M pc -2 (K km s -1 ) -1 when aco, in = 0.1 M pc -2 (K km s -1 ) -1 . 

In judging the biases of the actual solutions for our central solution pixels, we then select all of the aco, out mea- 
surements equal to aco,meas and find the mean and standard deviation of the aco, in values that gave those results. 
In Figure [l9| we show a plot of the mean and standard deviation of the aco, in values versus the aco,meas plus its 
uncertaintyior all galaxy centers in the sample. This figure shows that within the uncertainties of acomeas? none of 
the central pixel measurements show a bias. 
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Figure 18. Results of simulated data tests comparing the known input intrinsic scatter in DGR with the recovered minimum in the DGR 
scatter. The gray scale represents the number of trials falling in a given bin. We have uniformly sampled input intrinsic scatter values 
between 0.0 — 0.4 dex. The red line shows a one-to-one relationship. 
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Figure 19. This figure shows the results of simulated data tests focused on the galaxy centers, where high CO S/N and high CO/HI ratios 
could potentially bias the recovered aco- The x-axis shows the measured oico,meas and uncertainty for each of the 26 central points in our 
sample. The y-axis shows the mean and standard deviation of the known, input aico,in values from the simulations where we recovered 
«co,out equal to the value we have measured. This test shows that the aco values we measured in the galaxy centers are not biased 
within their uncertainties. 



In addition, a second check on these low central aco? a ^ least for NGC 6946, is shown in Figure 20 where we show the 
minimization curve for several techniques towards its central solution pixel. All techniques agree within their errors, 
highlighting the minimum at «co~ —0.4. This emphasizes that the low aco found in the center is not an artifact of 
the particular technique we are using, but can be recovered via minimization of scatter in the DGR, minimization of 
the correlation coefficient between CO/HI and DGR and minimization of the x 2 of the best-fit plane to CO, HI and 



Summary of Investigations into the Solution Technique 

To summarize, we have tested a number of variations on the basic technique of finding the most "uniform" DGR in a 
given region by stepping through a grid of aco values. We created simulated data with known aco and DGR intrinsic 
scatter matched to the range of S/N in our dataset. Via a Monte Carlo simulation, we found the robust mean of the 
logarithmic scatter to be the technique that spans the range of CO S/N and possible intrinsic DGR scatter with the 
highest accuracy and rate of success. Given knowledge of the CO S/N for a region and the measured minimum of the 
logarithmic scatter, we can use the simulated data tests to constrain how biased our aco measurement can be. We 
have also specifically checked that the central regions of galaxies, which only comprise 26 pixels out of the sample, are 
not biased due to their sometimes unusual CO/HI ratios. We found no significant bias for these solutions and show 
in addition that a variety of techniques can reproduce the low central aco solutions, not just the minimization of the 
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Log(«co) 

Figure 20. Solution for the central region of NGC 6946. The panels show different statistics we use to measure the uniformity of the 
DGR in a given solution pixel, all showing a clear minimum at olqo = 0.4Mq pc — 2 (K km s — 1 ) — 1 . A thorough discussion of the various 
statistics is presented in the Image Atlas. From top to bottom we show: the scatter in Log(DGR) measured using the robust biweight 
mean (LS.BIW); the fractional scatter DGR determined with the biweight mean (FS.BIW); the x 2 of the best-fit plane to Ico> ^Hl and 
Ed (PF); and the correlation coefficient between DGR and the CO/HI ratio (CHC). The dot-dash vertical line shows the solution for aco 
from the LS_BIW technique, which is what we adopt in the rest of the study. The agreement between different techniques (i.e. minimizing 
scatter in DGR, plane-fit and minimizing the correlation between DGR and CO/HI) illustrates the high confidence of this solution and 
shows that it is not an artifact of the solution technique. 

DGR scatter. 



Variations in R21 

In Figure |2l| we show the measured CO (2 — 1)/(1 — 0) line ratio in the solution pixels we have defined for each 
galaxy (the olack points) and in radial profiles (colored lines). The assumed R21 is shown with a horizontal solid, 
black line and a factor of 2 above and below that value are marked with black dashed lines. There are variations in 
but in regions with good CO S/N (the only regions in which we will achieve good solutions) these are generally 
less than a factor of two away from the R21 we have assumed. This deviation is within the uncertainties on the aco 
solutions, which typically have 0.2 dex uncertainty, at minimum. 
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Figure 21. The ratio of the CO J= (2 — 1) and (1 — 0) intensities measured from HERACLES and Nobeyama observations, respectively. 
The black points show all of the individual measurements from the solution pixels and the colored lines show radial profiles, binned by 
0.1r25- The value we have assumed for converting our (2 — 1) aco measurements to the more commonly used (1 — 0) scale is shown with 
a solid black line (R21 = 0.7), and the dashed lines show factors of 2 above and below that value. The variations of R21 are generally less 
than a factor of 2, which is less than the uncertainty on most of our aco measurements. Note that variations in R21 only affect the (1 — 0) 
acO) since we are directly measuring the conversion factor appropriate for the (2 — 1) line. 
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APPENDIX 

ATLAS OF GALAXIES AND RESULTS 

In the following Appendix we present the results for individual solution pixels in all of the galaxies we have considered. 
In the top panel of Figures [22] through Figure [46] we show the HI, CO and £d maps at matched resolution. The circles 
overlaid on these figures show the central position of each of the 37-point, hexagonal solution pixels in which we 
determine aco- In the middle panel of these same figures we show the aco values on the left in color and the 
associated uncertainties in the right in gray scale. Finally, in the bottom panels of the figures, we plot the measured 
aco values as a function of galactocentric radius (7*25). The color of the points in this plot reflects the uncertainties 
on the aco values as shown in the gray scale color table in the middle panel. 
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Figure 22. Results for NGC 0337. 
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Figure 23. Results for NGC 0925. 
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Figure 24. Results for NGC 2841. 
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Figure 25. Results for NGC 2976. 
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Figure 26. Results for NGC 3077. 
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Figure 27. Results for NGC 3184. 
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Figure 28. Results for NGC 3198. 



aco & DGR in Nearby Galaxies 

Log,o(N H ,) (cm- 8 ) Lo gl0 (l C0 ) (K km s' 1 ) 



Log 10 (S D ) (Mp pc- g ) 



45 



80.0 80.3 80.7 81.0 81.3 81.7 88.0 -1.0 -0.6 -0.8 0.8 0.7 1.1 1.5 -1.0 -0.8 -0.5 -0.8 0.0 0.8 0.5 



i ■ ■■ i ■ ■■ i ■ ■■ i ■ ■■ i 



i ■ ■■ i ■ ■■ i ■ ■■ i ■ ■■ i 



i ■ ■■ i ■ ■■ i ■ ■■ i ■ ■■ i 



48'0' 
46'0' 



48'0' 
40-0' 




11°36'0" - 



i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i 




i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i 




i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i 



10 h 44 m 80' 10' 0" 43 m 50* 40* 30* 



Log(q C o) 



Uncertainty Aa co (dex) 



Log (DGR) 



-1.0 -0.5 0.0 0.5 1.0 1.5 8.0 
■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ 



0.1 0.8 0.4 0.6 0.7 0.8 1.0 -3.5 -3.0 -8.5 -8.0 -1.5 -1.0 -0.5 



i ■ ■■ i ■ ■■ i ■ ■■ i ■ ■■ i 



i ■ ■■ i ■ ■■ i ■ ■■ i ■ ■■ i 



48'0" 
46'0" 
44'0" 
48'0" 
40'0" 
38'0" 
11°36'0" 




_i i i i i i i i i_ 





10 h 44 m 80" 10" 0" 43 m 50' 40' 



S 



2 - 



:mw 



o jo 



o 8 



o , o 

• v # %• 



o 
o 



o o 

■■■O-O 

O 



o : 

o o : 

co °; 

'0° G^ 



s ° o 



8 
o 



•• •■ 



^■■■•■■■O Q-a° O ^ean.^ 



(9 
o 



-1 - 



oo 



0.0 0.2 0.4 0.6 0.8 

Galactocentric Radius (r 86 ) 



1.0 



1.2 



Figure 29. Results for NGC 3351. 
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Figure 30. Results for NGC 3521. 
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Figure 31. Results for NGC 3627. 
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Figure 32. Results for NGC 3938. 
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Figure 33. Results for NGC 4236. 
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Figure 34. Results for NGC 4254. 
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Figure 35. Results for NGC 4321. 
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Figure 36. Results for NGC 4536. 
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Figure 37. Results for NGC 4569. 
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Figure 38. Results for NGC 4625. 
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Figure 39. Results for NGC 4631. 
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Figure 40. Results for NGC 4725. 
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Figure 41. Results for NGC 4736. 
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Figure 42. Results for NGC 5055. 
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Figure 43. Results for NGC 5457. 
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Figure 44. Results for NGC 5713. 
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Figure 45. Results for NGC 6946. 
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Figure 46. Results for NGC 7331. 



